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ABSTRACT 

An approximate analytic solution is developed for the problem 
of maximizing the range of an aircraft for a fixed end state. In 
general, this problem can not be solved analytically and is even 
very difficult to solve computationally. In this dissertation the 
problem is formulated as a singular perturbation problem and solved 
by means of matched inner and outer asymptotic expansions and the 
minimum principle of Pontryagin. 

Attention is focused on cruise at constant Mach number in 
the stratosphere, and on transition to and from cruise at constant 
Mach number. The state vector includes altitude, flight path angle, 
and mass. Normal acceleration and maneuvering drag effects are 
included. Lift is the control variable. Since Mach number is constant, 
thrust is constrained to be a function of state and control variables 
and is not itself a control variable. Specific fuel consumption becomes 
a linear function of power setting in the vicinity of cruise values. 

Cruise represents the outer solution. In cruise, altitude and 
flight path angle are essentially constant and only mass changes. In 
the inner solutions, corresponding to transitions between cruise and 
the specified initial and final conditions, mass is essentially constant 
and altitude and velocity vary. 

A solution is developed which is valid for cruise but which fails 
to satisfy the initial and final conditions. The cruise solution is 
shown to yield the Breguet range equation. By transforming the in- 
dependent variable near the initial and final conditions, we can seek 
solutions which are valid for the two inner solutions but not for cruise. 

The inner solutions can not be obtained without simplifying the 
state equations. However, to linearize them would completely elim- 
inate their dependence on altitude, as well as the dependence on alti- 
tude of any potential optimal control. The singular perturbation ap- 
proach overcomes this difficulty by allowing us to make a quadratic 
approximation to some of the state equations under certain circum- 
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stances. The resulting problem <s solved analytically, and the two 
inner solutions are matched to the outer solution. A modified Breguet 
range equation is developed which accounts for the changes in ramge 
due to starting from initial conditions not on a Breguet cruise and 
ending at final conditions not on a Bregeut cruise. 

The optimal control policy for transition is compared to several 
alternate control policies for supersonic cruise using the Boeing SST 
and for transonic cruise using the Boeing 707 and the McDonnell 
Douglas F-4. 
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CHAPTER I 


INTRODUCTION 


1. 1 Problem Statement 

This thesis is an attempt to develop an analytic solution to an 
airplane performance optimization problem that has proven very dif- 
ficult to solve by computational methods. It is the problem of maxi- 
mizing the cruise range of a gas-turbine powered aircraft. The 
problem is stated as a fixed end- state optimal control problem with 
a three dimensional state vector (altitude, flight path angle, and mass) 
in which the independent variable (range) is to be maximized. Atten- 
tion is focused on cruise at constant Mach number in the stratosphex-e. 
and on transition to and from cruise at constant Mach number. Simpli- 
fied computational solutions have indicated the existence of transonic 
and supersonic Mach number limiting, or constant Mach number cruise, 
for range optimal trajectories [l ] . By use of singular perturbation 
techniques and the minimum principle of Pontryagin, approximate 
analytic solutions are obtained as asymptotic expansions in three sep- 
arate segments of the problem: cruise, and transitions to and from 
cruise, and these segments are matched to form a single solution, 
valid from initial through final conditions. 

1. 2 Maximum Range in Cruise 

The problem of evaluating maximum range in cruise is of funda- 
mental importance to the design of any aircraft. It can not be separated 
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from the basic choice of mating a power-plant, with its thrust-producing 
and fuel -consuming characteristics, to an airframe with its lift and drag 
characteristics in cruise and its fuel-carrying capability. An airframe- 
power-plant combination will have one best altitude for cruise and one 
best velocity for cruise. At that speed and altitude the rate of consump- 
tion of fuel with respect to range is minimized. The aircraft must, of 
course, be able to maintain equilibrium flight at that altitude and velocity. 
The choice of the optimum altitude -velocity combination requires an 
iterative approach. At a trail point drag can be calculated. Referring 
to the engine characteristics (maximum thrust as a function of altitude 
and Mach number) one can determine if there is enough thrust available 
to balance the drag in equilibrium flight. If there is not, the trial point 
is not valid. For a valid trial point the fuel consumption, in pounds per 
second, corresponding to the required level of thrust is divided by velocity 
to give the parameter, pounds of fuel per mile, which is to be minimized. 
The minimum value will correspond to an altitude that is a compromise 
between the altitude for the most efficient unpowered flight by the aircraft 
(maximum lift-drag ratio) and that for the most efficient operation of the 
power-plant at constant velocity. 

Cruise velocity should be as large as possible but for gas-turbine 

powered aircraft it is limited by two considerations and hence will occur 

in one of two velocity regions, transonic or supersonic. For a transonic 

optimum the cruise velocity is limited by the beginning of the transonic 

drag rise. The rapid increase in drag associated with the transonic 

region translates to increased thrust required to balance the drag in 

cruise and to increased fuel consumption. For a supersonic optimum 

the cruise velocity is limited not by fuel consumption but by consideration 
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of the maximum temperature that the airframe can withstand [2 ]. 

Having determined a cruise condition, one can estimate the result- 
ing range capability by means of the Breguet range equation. This classic 
relation equates range to the product of a powered flight efficiency factor 
and the natural logarithm of the ratio of the initial mass of the aircraft 
to the final mass of the aircraft after its fuel budget is used up. The 
powered flight efficiency factor is a product of the lift- drag ratio (air- 
frame efficiency) and the ratio of cruise speed to specific fuel consump- 
tion (power-plant efficiency). This relation, of course, estimates range 
only for flight at the previously determined cruise condition and in no way 
accounts for flight to or from that condition. 

The Breguet range equation can be said to represent a one- 
dimensional approach to range capability estimation. To derive it one 
need consider only the state differential equation for mass, together with 
the equilibrium flight assumptions of lift equals weight and thrust equals 
drag. A correction factor can be derived to account for the fact that 
altitude does not remain constant in Breguet cruise but must slowly in- 
crease as fuel usage causes weight to decrease [ 3 ] . Only initial and 
final mass can be specified. Altitude is essentially a control variable, 
chosen at a particular value of mass to maximize the derivative of range 
with respect to fuel. 

Edelbaum has shown [ 4 ] that in the larger context of a maximum 
range cruise including initial transition to cruise and final transition from 
cruise, the Breguet cruise describes the optimal cruise portion for those 



problems In which range Is not so short as never to require a cruise 
portion. Edelbaum formulates a solution to the problem of range- optimal 
climb to cruise and descent from cru^e In terms of the energy- state 
method, another one- dimensional approach. Use of the energy- state 
method, described by various authors [1, 4 through 9], permits changes In 
velocity and thereby a complete solution from sea level to cruise. In the 
energy- state method normal acceleration is neglected (lift equals weight). 
As a result drag Is a non-linear function of altitude, velocity and mass. 

A recent study by Teren and Daniele [ 10 ] has expanded the state 
vector of the range- optimal cruise problem to two dimensions, altitude 
and mass, while neglecting normal acceleration and holding velocity con- 
stant. Thrust coefficient Is taken as the control variable. Again, lift 
equals weight and drag Is a non-linear function of altitude, mass and con- 
stant velocity. The problem Is formulated as a non-linear two- point 
boundary value problem and an approximate graphical method of solution 
Is presented. 

Kelley, Falco and Ball [ 11 ] used a four dimensional state vector 
(velocity, altitude, flight path angle and mass) In studying various airplane 
performance problems Including the maximum range problem. They were 
investigating the usefulness of the method of gradients in obtaining compu- 
tational solutions to these problems. For short range problems their op- 
timal result was a boost- glide or bang- bang solution. They reported that 
attempts at solving long range problems, which would Include a constant 
velocity cruise segment, were frustrated by convergence difficulties. 
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This thesis considers the range- optimal cruise problem with a 
three-dimensional state vector comprised of mass, altitude, and flight 
path angle, with lift becoming the control variable. Velocity remains 
constant, and so the problem is restricted to cruise and transitions to 
and from cruise at cruise velocity. Inclusion of normal acceleration 
(equation for flight path angle) means that the effect of lift, as well as 
altitude and velocity, on drag is included. Approximate analytic solu- 
tions are developed through the use of singular perturbation methods 
which, as will be shown in Chapter V, allow the drag force to be ex- 
pressed as a quadratic function of altitude, mass, and maneuvering 
lift. The equation for mass then becomes a quadratic function of altitude, 
flight path angle, mass and lift. The other two state equations are lin- 
earized and the resulting optimal control problem is solvable. 

1. 3 Singular Perturbation Problems 

A singular perturbation problem [ 12, 13, 14 ] can be described 
as a set of differential equations involving a small dimensionless parameter, 
say 6. The nature of the c- dependence is such that if c were to approach 
zero the order of the set of differential equations would be reduced. As a 
result the boundary conditions associated with the equations could not all 
be let simultaneously for a zero value of (. Viewed in another way one 
could say that the method of ordinary perturbations, involving the expansion 
of the dependent variables in power series in c. would fail to produce a 
solution that would be valid in the neighborhood of the boundary conditions. 
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Such problems are often solved by the method of matched asymptotic 
expansions, in which a stretching transformation applied to the independent 
variable in the neighborhood of the singularity transforms the problem to 
one that can be solved by ordinary perturbation methods in that vicinity. 
Then these solutions which are valid only in the neighborhood of the singu- 
larity can. by a choice of constants, be matched with those solutions that 
apply everywhere except in the neighborhood of the singularity to produce 
a single solution that will be valid throughout the region of interest of the 
problem. 

The range- optimal cruise problem seems well suited to formulation 
as a singular perturbation problem. It is convenient to think of the prob- 
lem as separable into a climb, a cruise and a descent. Certain variables, 
such ac altitude and flight path angle.undergo their greatest variations 
during climb and descent, but remain nearly constant during cruise. For 
mass the reverse is true, it being nearly constant in climb and descent, 
but varying most during cruise. Thus the problem is largely describable 
in terms of mass variation at nearly constant altitude and flight path angle 
except in "boundary layers" near initial and final time. We may think of 
mass as having its own characteristic time which is different from that 
of altitude and flight path angle. This characteristic of the problem makes 
it likely to be describable as a singular perturbation problem and offers 
the hope of yielding an approximate analytic solution that is uniformly valid 
over the entire time Interval of the problem. 



Interest in singular perturbation methods as applied to problems in 
aircraft dynamics begins with Ashley [ 15 ] . Drawing on an earlier work 
by Kevorkian [ 16 ] on reduced- order modelling, Ashley was able to dem- 
onstrate the separation of aircraft longitudinal dynamics into the short 
period and phugoid modes on the basis of the wide separation of their 
characteristics times. Kelley and Edelbaum [ 6] explored the idea of 
using singular perturbation methods to obtain a first order improvement 
to the energy- state solution to some optimal performance problems for 
airplanes. They thereby avoid the unrealistic instantaneous changes in 
altitude and velocity that occur in energy state solutions. Kelley also 
has suggested the use of singular perturbations in two-point boundary 
value problems and in reduced order modelling of aircraft performance 
problems [l7, 18, 19]. 

Kelley's objectives were to find reduced order approximations to 
certain airplane performance optimization problems that could still be 
related to the higher order computational solutions. These reduced order 
solutions could serve either to provide insights to improve the computa- 
tional solution or as good approximations in themselves to the higher 
order computational solutions. 

This thesis carries forward the ideas of Kelley and Edelbaum by 
setting up and solving the range optimal cruise problem as a singular 
perturbation problem. The perturbation parameter is developed naturally 
out of the parameters of the problem. The resulting solution is easily 
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related to lower order solutions and the nature of the solution gives in- 
dications of why the computational solutions are difficult to obtain. 

1. 4 Chapter Summary 

Chapter II presents and solves a problem similar in form to the 
state equations of the range optimal cruise problem. The solution of 
this problem demonstrates the techniques of solving a singular perturba- 
tion problem by means of stretching transformations applied to the in- 
dependent variable in the vicinity of singularities (boundary layers) and 
matched asymptotic expansions. 

Chapter III shows that the range optimal cruise problem can be 
expressed as a singular perturbation problem with singularities occurring 
at the initial and final state. 

In Chapter IV a solution is obtained to the problem of cruising 
flight which is valid everywhere except in the vicinity of the singularities. 
It is shown to be the Breguet solution. 

In Chapter V a solution is developed that is valid In the vicinity of 
general initial and final conditions but falls to be valid elsewhere. The 
solution is obtained by applying the minimum principle after expressing 
the problem as a linear optimal control problem with a quadratic cost. 

An optimal control is obtained and the optimal state trajectories are 
matched asymptotically to the cruise solution. A corrected Breguet 
equation is developed, accounting for fuel penalties (or bonuses) 
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associated with achieving Initial and final conditions. 

Chapter VI presents a cost comparison of the optimal trajectory 
with trajectories using various non- optimal controls In the near vicinity 
of cruise. It also studies the control as a sub- optimal control over large 
changes In altitude. Applications to several different aircraft are 
discussed. 

Chapter VII presents the conclusions and contributions of this thesis 
and suggests possible future work related to the thesis. 



CHAPTER II 


SINGULAR PERTURBATION PROBLEMS 


Singular perturbation problems and the techniques of solving them 
are most easily presented by formulating and solving a demonstration 
problem. The demonstration problem used here Is chosen for similarity 
to the state equations of the range optimization problem. It is adapted 
from O'Malley [ 13 ] . 


Consider the following set of differential equations in x, y, t and c 
X » y (2. 1) 

<y • -X - y (2. 2) 

where c is a small parameter. Initial values of x and y are specified 

x(t»0) « a (2.3) 

y(t«0) « b (2.4) 

It happens that this set of equations could be solved directly in terms 
of c to give the result 

+ r,t +rot 

(2.5) 


+ r.t +r«t 

X(t) ■ CjS + C 26 


+ r-t +r«t 

y(t) « TjCjS + TgCgC 


( 2 . 6 ) 


where 


r, • 




(2. 7) 
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i 

I 


i 


i 


I 


I 


i 

i 






■^ 11-^1 


(2.8) 


^ 2 * - 
^2- ri 


(2. 9) 


•r^a + b 
T2 - 


( 2 . 10 ) 


In cases where the equations can not be solved directly a useful 
technique is to assume that the dependent variables can be expanded in 
power series in c 


x(t, c) = XQ(t) + cXj(t) + c X 2 (t) 

y(t, e) » yQ(t) + cyj(t) + 


( 2 . U) 


( 2 . 12 ) 


where the and y. are functions of t only and are independent of c , 

At this point some terms will be defined as they relate to this thesis. 
More rigorous definitions of theso terms can be found in [ 14 ] . Consider 
the sequence 


2 3 

€, € , C 


(2. 13) 


As € -* 0 each term is small relative to the one preceding it. Such a 
sequence is called an asymptotic sequence. Our expansion for x(t, c) 

(and y(t, c)) is a series of functions of time weighted by successive terms 
of the asymptotic sequence in powers of f . If, as c “• 0, each term of the 
expansion is small relative to the one preceding it, the series can be called 
an asymptotic expansion. If in some domain of interest, D, the series is 
an asymptotic expansion for all values of t within D, then the asymptotic 
expansion is said to be uniformly valid within D, 


\ 





Now our solution technique calls for the asymptotic expansions for x 
and y to be substituted Into the original equations. The resulting ser;es 
expansion of the left hand side of each equation must equal the series ex- 
pansion of the right hand side of the equation. Since the x. and y. are In- 
dependent of c and since the expanded equations are valid for all small c, 
the coefficients of a given power of c on both sides of an equation must be 
equal. The zero order problem (coefficients of c^) Is solved first. The 
first order problem Is then solved In terms of the zero order problem, 
etc. This Is the approach of ordinary perturbations. For It to be suc- 
cessful the resulting solution must be uniformly valid In the time domain 
of Interest. 

Proceeding with our example, the expanded equations, to first 


order In c, are 

xq * yo *yi 
^y© * '<’‘0 yi' 

The zero order problem Is 

Xq » yg (2.16) 

Xq + yQ « 0 (2. 17) 

XQ(t«0) « a (2,18) 

yQ(t»0) » b (2,19) 
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Its solution Is 


Xq ■ ^36“^ (2.20) 

y^j • *' 3 ®** ( 2 . 21 ) 

The first order problem Is 

Xj » yj ( 2 . 22 ) 

Yq ■ •’^1 yi <2.23) 

Xj(t»0) « 0 (2.24) 

yj(t»0) • 0 (2.25) 

The result for y^ is now used in solving the first order problem 

Xi + y^ « -C3e‘* (2.26) 

Differentiating and substituting for x^ 

yj + yi » Cge'* (2.2V) 

yi » c^e”* *► Cjte"* (2.28) 

Xj ■ -c^e"* - Cge"^l+t) (2.29) 

Summarizing for x and y we have a solution to first-order in e or 
to two terms of the expansions 

X • C3®** ■ clc^e'* + Cge’Nl+t)] (2.30) 

y • -Cge"* + c[c^e'* Cgte’* J (2.31) 
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Now It It obvious that for c ■ 0 this solution cannot satisfy general 
initial conditions on both x and y. It can only satisfy one of them, and re> 
quires that the other be equal to it. 

a • -b (2.32) 

For any other initial conditions, one of the variables, say y, mue; 
make an instantaneous jump at t « O*** from its true initial condition to 
the value -x(0) and that jump is not necessarily of order c. That is, at 
X « 0 v'e must have y • so that we cannot say cy "• 0 as c 0. Borrow- 
ing a term from similar problems in fluid mechanics refer to this 
singular region as a boundary layer. Tl)e width of this boundary layer is 
of order c. Outside of this boundary layer, that Is, for 

t> c a 0 (2.33) 

the solution is valid. 

Such problems, for which ordinary perturbations fail, comprise a 
large class of singular perturbation problems. A direct way of identify- 
ing this kind of singular perturbation problem is by the fact that for < * 0 
the order of the system of first order differential equations is reduced. 
Instead of two differential eqtiations we have one algebraic and one dif- 
ferential equation. That means a reduction in the number of constants 
available for meeting boundary conditions and hence a failure to be able 
simultaneously to satisfy all boundary conditions. 

In order to analyze the region in which our solution fails to be valid 
it is useful to "stretch** the independent variable by the transformation 
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} ' ' 


r 

c 


(2. 34) 


with the differential relationship 


1 A . Am*) 

c dr dt 


(2.35) 


We can think of this device as allowini; us to view the problem on a faster 
time- scale as. for example, one would change the time- base on an oscil- 
lo^ t.ope to reveal an initial transient in what had appeared as a square 
wave. The transformed equations are 


A X » cy 

Or 


(2. 36) 


— y = -x-y 
dr 


(2.37) 


Again we expand dependent variables in powers of c and group 
corresponolng powers of c 


— X =0 
dT 0 


(2.38) 


Ay 

dT^O 


-*o ■ yo 


(2,39) 


Tr“> ■■ 


( 2 . 40 ) 




(2,41) 


f: i 
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The initial conditions are 


Xq(t * 0) » a (2.42) 

y^ir * 0) « b (2, 43) 

x^(t» 0) « yj(T* 0) « 0 I » 2. 3. ... (2.44) 

Notice that the transformed equations comprise a regular perturba- 
tion problem. In solving this problem it will now be possible to satisfy 
the boundary conditions. It is also important to notice that the transformed 
equations can be solved. If they cannot be solved, or if the untransformed 
equations cannot be solved, then the boundary layer transformation is of 
no value, since solutions of both problems are required. 


The zeroth order problem is solved directly as 


yo = 

The first order problem is solved as 



(2. 45) 

(2.46) 


(2. 47) 


Xj = -kje"^ - VCqT + kg (2.48) 

yj « -(kQ + kg) + kpT + kjTe"'*’ + kge"”^ (2.50) 
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Now it is possible to satisfy initial conditions to zeroth order in r. 
From the zeroth order solution we have 

kjj » a (2.51) 

kj = a + b (2. 52) 

The first order solution yields 

kg » a s b (2. 53) 

kg * 2a + b (2. 54) 

Now we have a solution to first order in c, or to two terms of the 
expansion, that is uniformly valid in the boundary layer. 

X = a + c[a(l - T - e"”^) + b(l - e~^) ] (2. 55) 

y = -a(l-e^) + be^ + c[a(-2 + t + 2e ^ + re ^) 

+ b(-l + re"’’' + e“^)] (2.56) 

The problem now is to reconcile this solution, valid in the boundary 
layer, with the previous solution, valid everywhere in the region except 
the boundary layer. It is common to call the boundary layer solution the 
inner solution. The other solution, which is valid everywhere in the region 
of interest except the boundary layer, is called the outer solution. In the 
problem of maximizing range in cruise there will be one outer solution 
(cruise) and two inner solutions. The first inner solution will describe 
the problem in the neighborhood of the initial conditions. The second 
inner solution will describe the problem in the neighborhood of the final 
conditions. 
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In our demonstration problem the Inner and outer soluti ms are com- 
bined t give a single uniformly valid solution by the technique of matched 
asymptotic expansions. The solutions are not matched at a point as one 
would match boundary conditions. Matching is based on the notion that the 
inner solution, valid in the boundary layer, and the outer solution, valid 
outside of the boundary layer, must both be valid in some overlap region. 

The inner solution is now extended to a form that it approaches 
beyond the boundary layer. First the independent variable is transformed 
to that of the outer solution 

T = -j (2.57) 

The solution is then expanded in powers of c. The resulting expansion is 
called the outer expansion of the inner solution. 

The outer solution is now extended to a form that is approaches as 
it approaches the boundary layer from some large value of t. First the 
independent variable is transformed to that of the inner solution 

t = €T (2.58) 

The solution is then expanded in powers of c. The resulting expansion is 
called the inner expansion of the outer solution. By suitable choice of the 
undetermined constants of the inner solution it will be possible to make 
the inner expansion of the outer solution identical, up to a certain order of 
f, to the outer expansion of the inner solution. 


First we evaluate the outer expansion of the inner solution. Using 
the transformation (2. 57) we have 


i 

! 


I 






I 
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(2. 59) 


« a(l - t) + <fa + b ] (1 - 
■ -a(l - t) + [a + b + t) 

-€(2a + b) (1 - (2. 60) 

These equations are expanded in powers of €, and in the limit of small c 
the exponential terms are vanishingly small. Using notation similar to 
that of O'Malley we write the outer expansion to two terms (order zero 
and one in <) of the inner solution to two terms as a function of t and € as 

[x(t, c)‘^l = a(l - t) + f(a + b) (2.61) 

[y(t. c)‘^| = -a(l - t) - c(2a + b) (2.62) 

Now we evaluate the inner expansion of the outer solution. Using 
the transformation (2. 58) we have 


o2 ^ 

Cgc"*'’’- cc^e’’’- ccgd + 

(2. 63) 

o2 . 

-Cge’'‘'‘+ cc^e*^‘’‘+ e^Cgte'^'’’ 

(2. 64) 


Expanding in powers of e we can write the inner expansion to first order 
in e of the outer solution to first order in c as a function of t and c 

[x(t, f)°^l * Cg - fCCjT + (2. 65) 

(y(T. f)°^l ® -Cg + c(CgT c^) (2.66) 

Transforming back to functions of t we have 
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[x(t. c)°^l = Cgd - t) - f(Cg + C 4 ) (2.67) 

[y(t. c)°^] » ‘Cgd- t) + <c^ (2.68) 

Now, the condition for matching is that 

[x(t. c)°^] » [x(t. c)^^l (2.69) 

[y(t. » [y(t, c)^^] (2.70) 

This is accomplished if we select the constants c^ and c^ as 

Cg a a (2. 71) 

c^ a -2a - b (2. 72) 


Now we can proceed to write a composite solution for x and y in 
terms of t valid to first order in e throughout the region of interest. To 
begin with, this solution will be the sum of the Inner and outer solutions. 
However, that implies doubly describing the variables in the overlap region 
a'here matching takes place. To remove this effect we subtract out the 
inner (or outer) expansion of the outer (or inner) solution. Finally, for 
our composite soultion we have 

x(t, = x(t, c)‘^ + x(t, c)°^ - [x(t, c)*^] (2.73) 

x(t)^^ = ad - t) + f[a + b](l - e'*^^) + ae’* 

- cad + t)e'* + c[2a + b] e"*^* - a(l - t) 

- c(a + b) (2. 74) 
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y(t)^^ « -a(l - t) + (a + b] e'^'^d + t) 

-c[ 2a + b ](i - - ae'* - c(2a + b]e** 

+ cate"* + a(l - t) + c(2a + b) (2, 75) 

These solutions simplify to 

x(t)^^^ * ae‘* + c[ae’*^* - ad + t)e‘*] (2. 76) 

y(t)'^^ » (a + b)d + t)e"*'* - ae“*(l - ct) 

- c[2a + b](e'* - e"*'*) (2.77) 

It is evident that these composite solutions satisfy the initial con- 
ditions exactly. For other values of t the error between these solutions 

2 

and the exact solutions given by (2. 5) and (2. 6) will be 0(c ). 
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CHAPTER III 


THE MAXIMUM RANGE PROBLEM AS A SINGULAR 
PERTURBATION PROBLEM 

In this Chapter it will be demonstrated that the state differential 
equations of the maximum range problem can be formulated in terms of 
a small dimensionless parameter, e. and that in the limit as c approaches 
zero, the order of the problem is reduced. Some assumptions that are 
used to simplify the equations are discussed and symbols are defined. 

The state differential equations for altitude (h), range (r), mass (m) 
and flight path angle (y) are, respectively 


ft » V siny 

(3.1) 

r = V cosy 

(3.2) 

ill . -JL(SFC) 

(3.3) 

8 


y s 3. /-i; — cosy) 

V V mcr / 

(3.4) 


Equation (3. 4) Incorporates the conventional assumption that the component 
of thrust (T) in the direction of lift (L) is degligible [l ]. 

True air speed (v) is assumed to be a constant. The constant speed 
cruise condition is predicted by range- optimal energy state solutions [I ]. 

It will occur either at the transonic drag rise or at the maximum supersonic 
Mach number. This assumption means that the problem is restricted to 
cruise and to constant speed transitions to and from cruise. It also means 
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that thrust Is no longer a control variable. Instead its value is constrained 
to be such that v remains constant 


T « W slny + D 

Aerodynamic drag. D. is described as follows 


D « CjjSq 


(3.5) 


(3. 6) 


of 


The wing area is S. The drag coefficient, C^, is assumed to be the sum 

a profile drag coefficient, or drag coefficient for zero lift, Cj^ , and a term 

o 

proportional to the square of the lift coefficient. The proportionality factor, 
K, is the coefficient of induced drag. Both Cj> and K are functions of Mach 
number. 

,2 




+ KC 


D » C-, Sq + 


KL 


L 

2 


D • Sq 


0 Sq 

+ KW2(l -I- u)2 
Sq 


(3. 7) 
(3.8) 
(3. 9) 


The dynamic pressure, q, is 

Q ■ 4- (3. 10) 

and if we restrict our problem to the stratosphere we have an isothermal 
atmosphere and two simplifications result: atmospheric density, p, becomes 
an exponential function of altitude 

-5(h - h.) 

p « p^e (3.11) 

and the speed of sound becomes a constant. Mach number is therefore a 
constant in view of our assumption of constant v, and and K are also 
constants. 
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The control variable, u, is defined as 


u > — - 1 

W 


(3. 12) 


Specific Fuel Consumption, SFC, is assumed to be descrlbable as a 
function of thrust coefficient, Cp, where 


T 

Sq 


(3. 13) 


It has been shown in [ 4 ] that SFC is a function of Mach number, power 
setting, and atmospheric temperature. Since our problem is restricted to 
constant Mach number flight in an isothermal atmosphere, SFC depends 
only on power setting, or Cp. The nature of this dependence is shown for 
typical transonic and supersonic cruising aircraft in Appendix A. In this 
problem we assume that in the vicinity of the cruise value of Cp we can 
express SFC as a linear function of Cp 


SFC « SFC^ +( - C„^) 

c dC„ F Fc 


(3. 14) 


We can also write SFC as 


SFC » Cq + CjCp 


(3. 15) 


atid, defining a constant, y, as 
c. 




(3. 16) 


we have 


SFC ■ c^jl + y 


1 


'oj 


(3. 17) 


24 




The constant, y, Is a measure of the slope of the curve of SFC as a function 
of Cp. It recurs throughout the rest of this thesis In connection with the 
description and derivation of the optimal solution. 


We now define the following dimensionless variables 



M . 

R . £ 
r* 


(3. 18) 
(3, 19) 


(3, 20} 
(3.21) 


where 8 is che scale height of the atmosphere and the asterisks denote 
reference values. For mass the reference value will be the Initial value, 
m^. It Is now possible to express weight as 


W » W.(l - M) 


(3.22) 


and drag as 


D » W. 



(1 - M)^(l + u)^ 


(3. 23) 


The reference value of q Is the value that minimizes the expression 
for drag In equilibrium flight, that Is, when T « D and L » W (u « 0), 
Differentiating (3. 8) with respect to q and solving for q we have 



(3.24) 
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Since L ■ W and v la constant in equilibrium flight, this equation defines 
the altitude for maximum lift-drag ratio. Also, since W Is a ratio of two 
dynamic pressures, we can express it as 


Q 


i, 

2 ‘b 




Q . ^-Kh-h*) 

Q . 


(3.25) 


(3.26) 

(3.27) 


Thus, from (3. 19) and (3. 26) we have tied the reference altitude to the 
reference weight. The reference value of altitude Is the altitude for the 
maximum lift-drag ratio attainable at the initial value of mass. At that 
altitude we Imve 

AH « 0 (3.28) 

Q - 1 (3.29) 


Now If we express the state equations in terms of the dimensionless 
variables AH, R, M and y, and convert from time to R as the independent 
variable by dividing by the equation 


dR 
dt ■ 


u 

—T cosy 
r* 


we have 


dAH 

dR 


Pr* tany 


(3. 30) 


(3.31) 
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OJI . 

dR \0v2/\cosy / 


(3. 32) 


^ . v: Ji)secy 

V \ Wj /\ ^0 Sq/ 


(3. 33) 


Since Cq has *.he dimensions of SFC, namely inverse seconds, the 
quantity (v/Cq) is a distance. We can therefore define r " as 


(3.34) 


Velocity will be on the order of 1,000 ft/sec, and typical values of Cq are 
about 0. 0005 /sec [ 1, 10, 20 j, so r'' will be on the order of 2 x 10® ft. 

In the stratosphere, P* ^ can be taken as 20,800 ft [21 |. Therefore, 
8r*, or , is dimensionless and its value is on the order of 100. We 
therefore choose its inverse as our perturbation parameter. 


(3. 35) 


It will be shown in Chapter IV that this parameter can be related to the 
cruise flight path angle. In fact, for y » 0 we have 


■ [HM] 


(3. 36) 


cruise 


Finally, in the equation for flight path angle we shall use the follow- 


ing definition 


a • — s* 

ev^ 


(3. 37) 
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This parameter Is much larger than c . In fact It can be related to Mach 
number, ^ , and hence Is on the order of one. 

2 1/2 

- tt * ^*'air 

where Is the ratio of the specific heat of air at constant pressure to 
that at constant volume and has a value of approximately 1. 4. 

The state equations can now be written as follows 


dAH 

dR 



(3.3P) 



(3.40) 


(3.41) 


It Is seen that in the limit as c approaches zero the differential 
equations for y and AH become algebraic equations defining y and u as 
zero, and arbitrary Initial conditions on y and AH could not be met. That 
is, for some non-zero initial value, y would have to go to zero in a zero 
Interval of range. It will be shown later that optimality considerations fix 
the constant value of AH which would also have to be achieved in a zero 
interval of range. Thus by demonstrating a dependence on c, a reduction 
in order in the limit as c approaches zero, and an inability to match given 
initial conditions in the limit as c approaches zero, we have demonstrated 
that the state differential equations of the maximum range problem can be 
formulated as a singular perturbation problem. 
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CHAPTER IV 


MAXIMUM RANGE CRUISE AND THE BREGUET RANGE EQUATION 


We now develop the outer, or cruise, solution. If in the state 
differential equations we assume that we can use a serieu approximation 
to two terms for trigonometric functions of y, we have 

tan y - y (4. 1) 


2 

cos y = 1 - ^ 


(4.2) 


Using these approximations and substituting fcr T, W and D from (3.5), 
(3.22) and (3.23), the state differential equations become 


dAH 1 
" 3 ^ 




(4.3) 




(4.4) 


dM 


{La-M,..yc^Q 




(4.5) 


The optimal control problem is to find the control, u, that 
transfers the state (AH, y, M), which is defined by the above three 
equations, from a given initial value, (AH^, y., M.) to a fixed final value 
(AHj, y^., Mj) while maximizing the final value of the independent variable. 


I 


I 

1 

{ 

\ 
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We can write the variational Hamilt<mian as 


.JK* ■*‘7 ^ y ^ ) ■*'7 ^ ~ 

^O (4.6 

Expanding the state and control variables in powers of < we have 


^ (AH„ * CAH,) • ^ ‘n * '"^2) 



(4. 7) 
(4. 8) 


It is already apparent that 
= 0 and 

«0 " I ^0^ = ° 


d 

HIT^O ' 


0 


(4.9) 

(4.10) 

(4.11) 


which implies that 

Uj = 0 (4.12) 

also. Therefore 


30 


TO <“o^ '“l> = -^] [>*>'(>* 

^ J /I ik/r \l 1 o— It _L ” » 1 _i_ Tjf 


0]} 


1> = Vi-p ^ [Wo^ -^j L‘ ’y V’ ns 

, .- , f/ 

+ C tyj(l-Mo)[l+2y(l+-..0 Ll'l 

VIq ® 

a-M )2 , (1-M >2 (1-M )2,, 

^0 


"0 

“^0 «i u u 

■i-rf[-^^L'^(>*^)>y-^(‘*^l 

0 ^0 S ^0 ^0 

(4.13) 

The costate differential equations 

fa . 1 - /_ v^\ ! 


y (l-M)y 

■ yc^ 


o 

I- y (i+lli^a+u)2)] + Xj^,y [d-M)y +yc:^ (q 
Q o 


o 


ja-M)^(i-Ki)^][i 
^ [d-M)y +yc;^ ( q+^(1-M)2( 


K (q*^,i-m)2(1^)=)](^)(1^)(i - ^)} 

o 'A 


(4.14) 


(i - ItW „^>2)r, ^ _J^ <1:™ 

(1-M)^ 

H ■ " '■■ 


•M)y 


„-k(q 

o 


i^(i^)2)]^XM[a-M)r*yc; 

<4..s) 
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(Hf) > 


+ y (l + (l^) (14^)2)] + Xj^[(l-M)y +yr^ (Q 


o 


We wish to find a control, u, such that 


(4.16) 


= 0 


(4.17) 


Performing the indicated partial differentiation of (4.6) we have 




+ y (l + (l+u)^)] + 2y [d-M)y + (q 

Q o 


- ^ a-M)^i-«.)2)] ^ 0^ I 

Q 


(4.18) 


Now expand the state, costate and control variables in these 
equations in powers of c. In doing so the following relations are used 


^0 * “o " 'H " ® 


(4.19) 


-AH -AHo-cAHj 

Q = Qq + «Qj = e ^ = e - Qq " ^^”l> 


(4.20) 


and, as a notational convenience 


1-Mq 


(4. 21) 


A 



The following expanded equations result 


" "c ^HO ^ {"^Hl " " ^MO 

+ 2y (1 + x^)]} +e {-Xj,2 ' +^0^2 > " ^M1 

+ 2y (1 +x2)1 - Xj^q7j <1 + +.V (1 +x^))] 

o 

2 

^ 2^^Qo+Xj^oMj(l+2y(l+3x2)] .Xj^oAHj(l-lVlo)l4yx^l} 
^^^o (4.22) 

" ^MO ^ Qo[a-x2)(l+y(l+x2)) 
o 

- 2yx^a+x2)]+c{x^, Q J(l-x^)(l4^'(l+x^)) 

^ o 

- 2yx (1+x )j ^MO % ^*^1 L 

o 

- (l+x^)(l-Hy(l+x^))-8yx'^]+Xj^Ov^^77T: Mj{2x)(l+y(2+6x^))} 

° (4.23) 

o 

+ c v'Cj. K (l+2y(l+x^)] + Xj^qV, [l+2yt€x^y) 

o 

^MO ^”l ll+2y+6x2y] 

O 

- ^MO (r^)<=‘’'>l>'^2y+6x2yl} 14. 24) 

O 0 

C‘V>l^’).l + €®(.X^;)o = c’WyO + 

+ Xj^o (l+2y(l+x^))]} (4.25) 
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Similarly we can express the expanded Hamiltonian to first order in c as 
•»'- • {-1 + + (1+x^) inydtx^)]} 

+ < {x„„y2 “So (“2 hti «o 


1 

I ‘‘ 


+ XMoyi<i-Mo)[i+2ya+x^)l + Xj^Q -/rp? QqAHj ((x^-ini+yid+x^)) 

O 

2„ . 2 


2yx^(l+x^)l - Xj^qMj (2xyC^) [l+2y (1+x^)]} 


(4.26) 


Now from (4. 25) we must have 


^yO ■ » 


(4.27) 


Also, in order to satisfy the expanded canonical equation for X we 

y 

must have 


^HO * ° 


(4.28) 


For these costates to remain zero over a non-zero interval of range, 
their derivatives must also be zero. Consider first Xjjq 


' 0 


(4.29 


TO ^HO " (d-x^Hl-tyd+x'')) - 2yx2(l+x2)) (4. 30) 


0 » 3yx^ + (l+2y) x^ - (l4y) 


(4.31) 


The positive real solution to this equation is 


X ■ + 


-l-2y +yi+16y+16y^ 
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(4.32) 


V 




^ r^^'-^.>r--^^--y if- s %-.t- '-« '-? 


It also will be useful later to express y as a function of x 

y * 


(1 - x^) 


(4.33) 


(3x^ - 1) (x^ + 1) 


Figure 4.1 is a plot of x vs y for physically realizeable values of y. 

Since y is the product of Cj^ and the slope of SFC vs Cp., y can never 

o 

be negative. Negative Cj^ is impossible in any case and a negative 

o 

d(SFC)/dCp would result in a "chattering" solution; the engine could 
be cycled on and off in such a way that its duty cycle would achieve 
minimum SFC. Therefore, x remains less than or equal to one, and 
other branches of this curve have no meaning for this problem. 


The fact that x remains less than or equal to one means that the 
optimal initial cruise altitude is below or at most equal to the optimal 
glide altitude, which is the altitude for maximum lift-drag ratio. This 
is seen from the definition of x when R is zero 


X = (1 - Mq-I/Qq. = 1/Qq. = e 


+AH 


Oi 


(4.34) 


or 

AHni = A«x -^ < X s 1 (4.35) 

0^ /T 

and from the definition of AH as the normalized altitude difference 
measured from the altitude for maximum lift-drag ratio. Therefore, 
the value of x given by Eq. (4.32) specifies an altitude at which the 
product of 1 /(SFC)q and lift-drag ratio has been maximized (to zero 
order in c) and this is lower than (possibly equal to) the altitude that 
corresponds to maximizing only the lift- drag ratio. 
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Using the optimal value of x, a constant, it is possible to solve 
the differential equation for Mq 


^ Mq = (1-Mq) (x + 1) (i-fy (i+x^)] 

0 

(4.36) 

= (1 - Mn)2xyUrTC 1 

" o Ox'^-l) 

(4.37) 

-2xyrpr(l+x^) R/ (3x^-1) 


1 - Mq = ce ° 

(4.38) 

and the constant, c. is 1 so that Mq is zero when R is zero. And now. 

since 


-AHq 

1 - Mq = xQq = xe " 

(4.39) 

we have 


AHrt * 2x K R + X 

° ^o (3x'^-l) 

(4.40) 

which implies that 


^ °o (3x'^-l) 

(4. 41) 

which is a constant, implying that 


1 . 2 
«2 “ ’7^1 

(4.42) 

Notice that initial values of Vq and AHq are fixed. 

and cannot be matched 

to arbitrary initial conditions. Alternatively it could be said that from 
arbitrary initial conditions Vq and AHq must move to their optimal values 
in a zero interval of range, thus demonstrating the singular nature of the 


problem in the vicinity of initial conditions. 


37 


Finally, the equation for Mq can be written as 


1 ■ Mq = e 


•y^R 


(4.43) 


Another way of defining y^ is as follows 

>1 ■ ^ (* ■•■ 7) + y (1 


(4.44) 


.DRAGq^ , SFCq^ 

" LtifTq'J L~c^J 


(4.45) 


This will be useful later in developing the Breguet equation. In this 
form and (SFC)q refer to the zero order problem, in which lift 

equals weight and thrust equals drag. Since the cruise flight path angle 


is cy^ to first order in € we see that € can be expressed as 
(sFC^) [y ^^^0^ BREGUET 


(4.46) 


If y=0, this expression is further simplified by the fact that Cq=SFCq. 


Concluding now with the zero-order problem we have from the 


condition that - 0 that 


- 1 + Xmo + j) [1 + y (1 + x^] = 0 


(4.47) 


- 1 ^Mo n 

/ 1 \ 

%0 = ip ^ 

This result is consistent with the canonical equation for Xj^q which 


states that 


m ^Mo 


= Xj^Q 2x ^/CJ 3 K [l+2y (l+x^)I = yj X 


■MO 


(4.50) 
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Now using the optimal constant value of x and resulting expression 
for Mq it is possible to solve the differential equation for 

^ Mj = Vj (1 -Mq) [l+2y (1+x^)] + AHj Qq [{x^-D (1+y (1+x^)) 

+ 2yx^ {l+x^)l + l-2x^ (l+2y (l+x^))) (4. 51) 

From the optimality condition on x the coefficient of is identically 
zero. Using the definitions of and 1 "Mq the equation becomes 

2 

°M. = i e ^ - y.M. (4.52) 

^o 


+ y R 

An Integrating factor is e and we have 


^ 2x k 

o 


dR + c 


-y.R - , y, R , ^ 

o 

The constant of integration is not necessarily zero. It represents a 
first order initial value of mass in cruise. 


At this point is is possible to develop the Breguet range equation. 
However, first we shall investigate the first order necessary conditions 
for optimality and determine AHj, y 2 and Ug. The condition that the 
derivative of X^q must be zero requires that 


- Xfji - Xj^o (l-Mg) (l+2y (l+x-^)] = 0 


2xyCjj K 


(4.56) 
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This is a constant and hence its derivative must be zero 


o 

■ ^MO^l % <l-*y a+x^))-8yx'*l 


+ Xj^O 2x,yu^ Mj (l+2y (l+Sx^)] 


(4.57) 


Here the coefficient of Xj^j is equal to zero at the optimal value of x. 
0 ^ i + AH, Hl+x^) - ya+2«^+9x'*)l 




2 2 
+ 2x^ (142y (l+3x^)) 


(4.58) 


0 = -4yjX^ (1-x^) - 2x^ AHj (-3x^ + 6x^ + 1) 

o 

+ (2x^ v'Cjj K Mjjj(0)+ X (-3x^ + 6x^ + 1) (4. 59) 


AH, = 




^ 2xyCjj k K (l+6x''^-3x^) 


2x y,(l-x^) 


(4.60) 


This Implies that 


^2 


2x7^^ 


(4.61) 


Again this is a constant. From comparing second power terms in the 
expanded differential equation for y, we must have 


Hi^y2 » «3''V2 


(4.62) 
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'1 

u, = i (4.63) 

^ 2xyu^ 

o 

It is now possible to determine from its equation 

d 2 / D^-i \n 

HlT^Ml ^MO <^‘*'2y+6x y) [_yj+ 2x K 6^) 

o 0 

- %R 

An integrating factor is e , and we have 

+ %R 9 ^ 

^M1 " ^ ^ (l+2y+6xV) [1 

2xVCd R 2x 


(-— ) 

' /7^ TT/ 


(l+6x^ -3x^) 




^M1 ® 


+ VjR 


1 -L r> 


] 


(4. 65) 


(4.66) 


Ox*”-!) 

Now the first order term in the Hamiltonian becomes 

o 

+ (i^^fl+XMod-Moiyc^ <> 

o 

+ y (1+x^)) + 2xy (1+x^)] - Xj^qMj 2x [l+2y (l+x^)l (4. 67) 

o 

The coefficient of AHj is zero at the optimal value of x, and the 
remaining coefficients are more conveniently expressed in terms of 

Substituting for Xjij, Xj^j. 1 -Mq, Xj^q. and from (4.56), (4.66), 
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4 






I 


(4.43), (4.49), and (4. 54) we have 


»'r- 




2x.yCp k ^ 2xyCp K"^ 2xyCj^ k 2xyCp k 


Mjjj(O) (4.69) 


•^1 ' 


(4. 70) 


and for = 0 we have 




(4. 71) 


This constant originated in (4.66) so we now have for \j 


[- 1 ^— +-J^1 

2 x yCp K 


(4. 72) 


Finally, can be determined from the condition that = 0 as 


follows 


ax 1 + Xj^o ( 1 -Mq) 2xyr^ (l+2y (l+x^)) = 0 
^ o 

and, substituting for Xj^jq, Mq, and we have 


(4. 73) 


(4. 74) 


We have now satisfied the necessary conditions for optimality to 
first order in c. Figure 4.2 is a diagram of the sequences that led to 
the zero order optimal solution. Figure 4. 3 is the s'me lor the first 
order optimal solution. An arrowhead from one box wO another indicates 
that the information in the first box leads to the conclusion in the second. 
A summation of two arrowheads indicates that two information sources 
are necessary to draw the indicated conclusion. It doe 7 not indicate a 
summation of equations. 


S 




Figure 4. 2. Zero Order Optimal Outer Solution Sequence 
































We return now to the expression for M and develop the Breguet 


range equation. 


M, = (Mrt + cM,)^ = 1 - e ‘ +<e ‘i — — 

r ^1 "I 

1 - Mf = e ^ M 1 - f 1 — 

^ ^ 2x 

° 2 

r ^1 1 

i** (l-MJ = - y.Rf + 1 - € — — 

‘ 2xyc^-‘ 


(1-Mf) = - yjRj. - c 


Rf = - l» (1-Mf) 


2x yCp K 
*■ o 

1 


2xyr^' 


■] 


<y 

R = - (1-M.) (^1) ( 1 - c i ) 

f ^ 2yiJTr^' 


Now using the following relations 


. tr, - mps .m.v 


and 


1 _ ^0 


yj draCq yiT- 


and 


■ «f (^) 


we have 


(4.75) 
(4. 76) 

(4.77) 

(4.78) 

(4.79) 

(4.80) 

(4. 81) 
(4.82) 
(4. 83) 
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“ Ksirr „„ ^ 


0 max 


2x .'rTTK'- 


(4.84) 


D 


The first brackets contain the usual Breguet range expression. 
The optimal value of x maximizes the product of v/(SI'C)q and (L/D)q 
and hence maximizes range at constant velocity and fixed initial and 
final values of mass. 


The second brackets contain a first order correction factor 
resulting from the fact that flight path angle is not zero but a small 
position quantity and therefore thrust equals not only drag but also a 
component of weight in the thrust direction. This extra thrust required 
for climbing results in a smaller final range, but the difference is of 
order c. 


Rutowski [2] derived an expression for this range correction factor 
from consideration of the increase in potential energy due to climbing 
during Breguet cruise. He considered SFC to be constant in his develop- 
ment and did not attempt a mathematical maximization of the resulting 
range expression. If the above expression for maximum range were 
derived based on u constant value of SFC, then the resulting maximizing 

value of X would be unity and would be 2 k. As a result, the 

o 

correction factor appearing above would become (l-c) which agrees 


exactly with Rutowski's 


factor of 




Teren and Daniele [10] have analyzed the maximum cruise range 

problem using only mass and altitude as state variables and using thrust 

coefficient as the control variable. They derive the following equation to 

define the optimum value of C-p 
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^r— B 


(SFO(Cp) 


(SFCHCj,,) 


2 i(SFC)(Cp)l 2^v 


(•4.85) 


If SFC is taken as varying linearly with Cj, in the neighborhood of a 
constant operating value, this equation is a cubic in ^ and c: 

r\ 

2yC 


r Cp 

o 


o o o o 


(4. 86) 


For c = 0 the order is reduced by one and the positive real root is 
easily found as 

(Cp)Q = Cj^ (1+x^) (4.87) 

o 

For ( 0 a first improvement to this root is found by assuming it to 

be of the form 


Cf = (Cp.)o(l + 6) 


(4. 88) 


This is substituted into the cubic and solved for 6, retaining terms of 
0(c) and eliminating terms that comprise the second order equation for 
(Cp)Q. The result is 


Cf = (Cj, ) (Kx 2) [_1 f C ^ 


(x2.1)2 

T 


(4.89) 


o (Sx**-!) a+ex^-Sx’)-* 

The identical result is achieved if (T/Sq) is expanded using the 
optimal values of Mq, Mj, AHj, and developed in this chapter. 


Teren and Daniele do not develop an expression for maximum 
range in cruise. However this is easily done using their differential 
equations for mass and range together with the assumptions of constant 
lift coefficient and lift equals weight. The result is 
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( 4 . 90 ) 


where SFC and Cj^ are determined by the optimal value of Cp. 

Now if SFC is a linear function of Cp and if the above expanded 
form of the optimal value of Cp is used, together with the resulting 
expanded forms of the optimal C^ and SFC, this expression for 
maximum range will be identical to the one derived in this chapter. 


CHAPTER V 


MAXIMUM RANGE TRANSITIONS IN THE BOUNDARY LAVERS 


5. 1 Introduction 

In this Chapter the two inner, or boundary layer, solutions are 
derived, and the results are matched with the outer or cruise solution. 
The optimal control is examined and shown to produce a damped oscilla- 
tory transition to and from cruise. The over-damped, or pure exponen- 
tial, case is examined separately. The short-range problem, for which 
no cruise segment is required, is also examined. Finally, a modified 
Breguet range equation is derived which includes changes In range due to 
transitions between cruise and initial and final conditions that are not on 
a Breguet cruise. 

5. 2 The Problem in the Boundary Layer 

The boundary layer problem is described by stretching the indepen- 
dent variable in the state equations by the transformations 



in the left side boundary (vicinity of R = 0) and 


R 


Rf - R 


R 


(5.2) 


in the right side boundary (vicinity of R = R^). In the left side boundary 
layer ^ f and the equations become 
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1 




dAH . 

aiq: 


■L^ O 


o 


The cost to be minimized becomes 


R^/c 


- i ‘ <‘”l 


and the variational Hamiltonian becomes 

2 

+X a(u+^)+Xj^y + €\j^ {|jl-M)y + 7Cp k[Q 
' o 


(1-Mr 


(l+u)^]] [i 


(l-M)y ^ 

• T 


o 


* 11 ^ , 1 ^) 2 ]] ( 1 +^)} 


(5.3) 


(5.4) 


+ y [i + Ik^ (1 -Ki)2 J] (1 + ^)| (5.5) 


(5.6) 


( 5 . 7) 


Unlike the situation in cruise, a zero order analytic solution to the 
boundary layer problem can not be found unless some further simplifying 
assumptions are made about the state and control variables. Typically in 
such cases one might consider a linear expansion of the boundary layer 
state t^K^uations in the vicinity of cruise. There u, y and M will be small. 
Let ^ be defined relative to the initial optimal cruise altitude so that 
also remains small. 
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= Q/Qq. = Qx 


’ i /K 




AH = AH - £>/x 


The linear differential equations are 


3 ^ ==“ 


^ • y 


^ =€{v^C; 71< (X +i) (l+y (l+x^)] +y(l + 2y (1 + x^)] 


+ ^H K {2x^y (^ + x) + (x - ^H1 + y (1 + x^))) 


■» 2 K X (u-M) U + 2y (1 + x^))} 


The equation for M is derived by expressing the exponential form 

A A 

of Q as a Taylor series in AH to first order, expanding, and retaining 
only linear terms. It is then greatly simplified by making use of two 
algebraic identities. First, from the definition of the optimal value of 
X (Eq. 4.31, 4.32) 


2x^y (X + -) + (X - i) (1 + y (1 + x^)] = 0 


(5.14) 


1. 


Second, the first-order value of the cruise flight path angle is defined by 
the two equivalent expressions (Eq. 4.44, 4.50) 


^Ib =^/C^(x +i)[l+ya+x^)J = 2yC^x[l+2y(l+x^)l (5.15) 


o o 

The differential equation for M then reduces to 


^ (l + u - M + "I ) 

o 

Now the Hamiltonian can be written as 


(5.16) 


.If - - i aX.u + X,.y + cX«,y,.. (l + u - M + ) (5.17) 

y 2xyc^^ 

o 

A 

It can be seen at this point that AH does not appear at all in the 
state equations or the Hamiltonian. An optimal control, if one could be 

A 

found, would be Independent of AH. Thus linearization of the state 
equations fails to yield a meaningful solution because it fails to yield 
a meaningful mass equation. The reason for the failure is that the drag 
term, which appears in the mass equation, has no linear dependence on 
altitude in the vicinity of cruise altitude. 


As a second attempt to simplify the boundary layer state equations, 
and thereby to be able to develop an analytic solution, consider a quad- 
ratic expansion of the equations in the vicinity of cruise. The differential 

equation for AH remains linear (5.4). The equation for ynow includes a 
2 

term in y (5.3). The equation for M becomes 


52 


! 


HM r 1 ^1^2 A 1 A9 

€ \ [^(l-M)y+,/c;7K [i(l-AH+j^H'‘)+x(l + 4*II+^4FD(l 


M)^(l+u)^]][l (1 +LH +juH^)|_(l-M) y+%T^ ■ 


/Cj, K 
o 


+ j**H^)+x(l + uH+^AlV)(l-M)^(l+u)“]]] (1+^)} ^5.18) 

A A 

with Q expressed as a Taylor series in AH to second order. Pei forming 
the indicated multiplications and neglecting terms of higher order tnan 
second results in the following equation 

^ = c{yc;p<(x+i)[l+y(l+x2)] + y[l + 2y(l+ x^)) 

L o 

+ AH,/C^ (2x^y(x +i) +(x -i)[l +y(l +x^)])+(u-M) 2xv'cyjKH+2y(l+x^)] 




v'C„ K 


o ,._1 


o 


X +^)[1 +y (1 + X^)]) + yAH4x^y 


+ AH^ ,€^(4x^y+j (x+i)[l+y(l +x^)]) + yu4x^y 


+ (M^ +u^) X /Cjj K [1 +2y (1 + x^) + 4yx^] - y M [1 +4yx^ +2y (1 +x^)] 
o 

+ (AHu- AHM) 2x yc,. K[l+2y(l+x^)+4x^y] 


- Mu4x yCp K [1 +2y (1 + x^) +2x^y)| 


(5,19) 


These equations (5.3, 5.4, 5.19) are solvable as the state equations 
of a linear-quadratic singular perturbation optimal control problem if we 

y2 

can neglect the term in comparison to u in Kq. (5.3). in the cruise 
problem this term had no effect on the solution through first order in €. 


M 


To show qualitatively that this term is negligible in the boundary layer, 
consider the effect of the transformation 

U*u+^. (5.20) 


on the state equations. The equation for y (5.3) is linearized and the 

A 

equation for AH (5. 4) is unchanged. The effect on the equation for M 

(5.19) can be shown to be merely the replacement of u by U and the 

2 

inclusion of an additive term in the coefficient of y . This additive 
term can be shown to be so small in comparison to the principal term 
in that coefficient that it is safely neglected. But if that term is neglected, 
the transformed equation set can not be distinguished from the set that 
results from merely neglecting the term -hp in comparison to u in the 
original equation for y 


Of course, the entire solution of this problem could be based on 
this transformation and no appreciable difference would occur in the 
result. For our present purposes we merely cite the potential of this 
transformation and proceed to Imearize the equati''n for yby neglecting 

v2 

in comparison to u. The equation for M is simplified by the use of 
(5.14) and (5.15). Finally, the state equations of the linear-quadratic 
singular perturbation optimal control problem and the variational 
Hamiltonian are 



au 


(5.21) 


dAH 


y 


'5.22) 
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I 


dM 

-3TT 


i . ..... y- .) 4 ^g_) 

L ^ 2xv^Cj7K^ vC,j 


+ (^4yx^ ^--.- ^ J . L . , j^-yMH-2Xs Cj^ k(aH|u-M) +iL-lM^lAIL)J 
2x.^/Cj^ K o 


+ y (u+AH)4x^y-4xv'C,. K Mu ( — I- -- + 2x^y)} (5.23) 

^ o 2Xv/C’,^ K 


5x^/C j^ K 
o 


- < +ax. u +X„y 4 cX {r,„ ^ 4 „ . M 4 


2k , Cjj K' 
‘ o 


y 

+ +y(u+AH)4x^y - 4xv^Cjj R Mu ^2x^y + 






2x ..r'lj 

o 


« y 

-t ^4yx‘^ + — . - ^l^-yM + 2x /Cj^ K (ah [u-M] + 

2Xv'^Cj^ K o 


u^ +M^ +AH*‘^ 

2 ; ■; 


(5.24) 


The factor of c in the equation for M makes this a solvable problem. 
The c means that the equation for M is not in the Hamiltonian to zero 
order in € and, as a result, the zero order value of the costate for mass, 
^MO' constant. Thus the zero order problem, for which mass is 

A 

constant, has linear differential equations for y and AH and a quadratic 
cost functional which represents an "out of the loop" equation for M 
weighted by a constant, X.j^q> The first order problem is of course 
linear with coefficients depending on the zero order problem. 


We now proceed to develop the costate diffei ential equation 
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dX 


^ - c . y (y,, . - 0 =) . (u . aA) 4x^y 


2xyCjj K 


yCjP< 


- M 


(4x y 


lb 


2xyCjj K 


)} 


dX 

HTT 


~= -«XM{y4x2y+ (yi 5 + 8 x^yvc;p<) (u-M + ah)} 
L ‘ o 


(5.25) 


(5,26) 


HIC ' ‘ {■ ’'lb ■ ^ ’'lb ® ^ * y) 


+ ^ 4 x^y+ y + 2 k, fC^ (M - AH)]} (5.27) 


2 x yCjj K 


The condition that the optimal control, u, minimizes the 
Hamiltonian is expressed as follows 

•^u = ®^y ^ L’'lb ’'o • (^^Ib ■■ ^ 

^ o 

+ (yih + 8x^y k) ( AH + u)] = 0 (5.28) 

o 

Now if we expand the state, costate and control variables in pc^^er 
series in €, the differential equations of the state and costate variables 
become 


31^ <V0 + f Vj) = a (Uq + < u^) (5. 29) 

^ (AK^, + € AHj) = yo + < Vi <5- 30) 

1 ^ 
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+ ^4 x^YvCj 3 R +^)[uq +AhJ + 4 xV Vq (^0 
o 

3 ^(x ^ + a ^ c\ ^2 ) = - x„o - « (xhi + Xj,o Vo (no * J^) 

o 


o 


+ 4xV ,u„ 4 4ft„)]} V €2{- x„2 - X„, 

o 

+ Vo(nb +-p^) +4x^y (Uq + AHq)] - Xj^o[yi 

v^D ^ 
o 

+-^SL.y 4x^y (u, + AH, ) - M, (4x^v + — )!} 

jr^' ^ 1 1 IX 2 xvF^ 

O o 

(5.32) 


<^h , -""^111 >'o ■'(^b -^^^o)]} 




• ^0 ■" (nb ■" k)(uq + ahq)] 

o 

^Mo['**'“v V, + (v,b + 8x^ vC^) (v “i + ^Hl)]} 


(5.33) 


<Xmo ^'^ 1* ■ ■ ‘^MO {- Vlb ' '’’o * Vib) 


(5.34) 


Expanding the equation = 0 gi es 


\ 



+ < {*^y 1 ^Mot^lb '‘’ ^* y ^0 '*’ (^Ib '*’ y'^Dj^)(^^0 “o)]} 

+ «^{aXy2 ^^Mlkb >0 ^(^Ib +8 x^/T^)(aHq +Uq)] 

^Mo[^*^y ^1 ■ (^^Ib ■" 8x\ve^) Mj + (xib + +Uj)]} 

° (5.35) 

Similarly, the expanded Hamiltonian is 

•»’- Ho“o -"^HO^o} {- ' -^'^Hon *^Hl’’0> 


o 




;)^C:j^K x(Uq+AH(,)^]} {a(Xyo^2 


W4x2y+— 

2xyU^'’ 

o 

* iXho^'o *^H2’'o> *%i[nb(* +" 


fl+U„+ ' 

i.y2(2^._JSL.') 

° 2xyCj. K^ 
o 

0 


+ roibo + ^"o’ ^*^5' ■" * '"o ■" 

2x.^ p K o 
o 

" )^Mo[’'lb(’-l - “l (’'lb '’'“b 

* ■'■ ^’'ib) 

■l («^ + ^7^=^)L‘ VoM,*2xVCp^K (aHjUj 
+ Afijbj - AH„Mj + u„Uj + AH„AH,)]} 


(5.36) 
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Because is a constant we can now see that the costate equations 
are linear and the Hamlltonoan is quadratic and we should now be able to 
develop and analytic solution. 


5.3 The Optimal Control 


The required optimal control, u, must minimize the Hamiltonian, 


that is 


0 


(5.37) 


Considering first the zero-order part of we must have 




(5.38) 


and since the Hamiltonian must itself be zero along the optimal trajectory 
we must also have, from = 0 


^HO " ° 


(5.39) 


Now the optimal Uq can be determined from the first order part 


of .yf, = 0 




(5,40) 


* (^Ib (o„ + AH„)] = 0 (5. 41) 

o 

An expression for the optimal value of u^ results from considering 
the second order part of = 0 




(5.42) 





(Jlf„)j ■ »Xy2 + Xj„[yi^+4x^yo +(y^, + 8x^vP^)(u(, + AH„)] 

o o 

(5.43) 

'.Ve can now solve for the optimal and u^. It will be shown that the 
optimal Uq satisfies a homogeneous fourth order linear differential 
equation with constant coefficients. Furthermore, the coefficients of the 
first and third derivatives are zero, indicating that the roots are symmet- 
rically located in the complex plane with respect to both the real and the 
imaginary axes. Next the optimal Uj will be shown to satisfy a non-homo- 
geneous equation, the homogeneous part of which is identical to the 
equation for the optimal Uq. The non-homogeneous part, or forcing part, 
will be a function of Mj which is itself a function of the optimal trajectory 
of the zero-order problem. 

First, since (dJ(''/5u)j - 0 over the entire optimal trajectory, it 
follows that the derivative of this partial derivative, with respect to R, 
the independent variable, will also be zero over the entire trajectory. 

Such derivatives will be indicated by primes. Similarly, all higher order 
derivatives will be zero. Proceeding, then, to tak ccessive derivatives 
of (d.3lfl'au)j, the differential equation for the optimal Uq is derived. At 
each step the derivative of a state or costate variable is replaced by its 
defining canonical differential equation. 

y 0 ■^(^ib 0 o)] 

L- ' O 
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(5.45) 


0 ■ .[- + >0 (y,^ *^=) * “■‘V AHoJ] 


•2xyc^ 




* ^MO (’'lb ■*■ ’/^'d *^)(’’o ■*' “ o)] 


dS, 




(’o *"'o) 

o 

° * ^^Mo(yib ./Cp 

o 

* ^Mo (®*^y * nb)(““o ■" “'o) 


(5.46) 


(5.47) 


No further substitution for variables is required and the next two 
derivatives are taken as a single operation 

^'ib 

2 


IV . 

“o ■"'*0 


V.K+-^ 

Tcpr 


2a - a‘ 


Vjb+8* 


D. 


+ a Uq = 0 


(5.48) 


Using p as a derivative operator and defining Z as 


v,b*^ 


Z = 




x/C^ 


(5.49) 


Vjb+8x yVCjj^K 


the differential equation for Uq may be written as 

‘o 


[p^ + p2 (2a - (2Z/al^) + a^J u. = 0 


(5.50) 


or in the equivalent factored form 
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(p^ + 2Z7a p + a) (p^ - 2ZyiT p + a) Uq = 0 


(5.51) 


The roots of this equation depend on a, which is a function of v, 

and Z, which is a function of v, x, y, and K. At a given Mach 

o 

number and K are constants. Since x is a function of y we can say 

o /d(SFC)\ 

that the roots depend on velocity and y, that is, (— . If at 

'“F o 

a given velocity d(SFC)/dCp is allowed to vary between zero (constant 
SFC) and + •, the root locations for a particular airplane will vary 
as shown in Fig. 5.1. It is obvious that the effect of increasing the slope, 
d(SFC)/dCp, away from zero is to increase the damping ratio. The 
natural frequency on the oscillatory brench remains constant. 

A locus of roots as velocity varies is somewhat more difficult to 

obtain because of the fact that y as well as Cp and K would have to be 

o 

described as a function of Mach number. However, if y is zero, as is 
sometimes necessary to assume, this difficulty is greatly diminished. 
Such a locus of roots as velocity varies and y remains zero is shown for 
a particular airplane in Fig. 5.2. In this figure it is seen that the effect 
of increasing velocity is to reduce the damping and the natural frequency 
of the oscillatory branch. Below a velocity of 409 ft/ sec the roots are 
non- os c illato ry . 

Returning now to Fig. 5.1 and the equation for Uq we can say that 
on the oscillatory branch we have for natural frequency and damping 
ratio, respectively, 

(5.52) 


C = z 
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(5.53) 







■-4 


When the roots of the equation are all real they are 

p = ±v^[z - l] (5.54) 

The damping is principally attributable to the slope of the curve of 

SFC vs Cp for the engine. From (5.51) we see that if Z were zero, the 

dynamic modes represented by that equation would have no damping. The 

damping comes from the term a^ + 2xy/ K^. Each of these 

2 ® 

terms originates in a term in the cost functional. Since y^ is propor- 
tional to dAH/dRj^, they represent a penalty on altitude rate excusions. 

The first of these terms can be traced to the effect of the cosine of 

the flight path angle on range. The second comes from a combination of 

the restriction of velocity to be constant and of the modeling of SFC as a 

2 

linear function of Cp. Mass rate is therefore proportional to T , and 

2 

since thrust must have a component equal to Wy, the second y term 
enters the cost functional. The second contribution to damping is the 
larger one. 

The equation for the optimal value of Uq can be solved to express Uq 

A 

as a function of Rj^. Then, in turn, the equations for y^ and AHq can be 
solved and the entire zero-order state and control are known as functions 
of Rj^. Finally, it is seen that Mj, Xj^j, X^j, and Xj^j are functions of the 
zero order state and control and hence they, too, can be expressed as 
functions of Rj^. 

Returning now to the expression for the optimal value of Uj, we can 

say 
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where 


. i 


("l) ' ■ ^Mlhb ■"■'*^3' ’’O (^b yCp^K)(uj + AH„)] 

o 

or using (a.X'T'au)^ = 0, 

^1 (^l) (x^) ^MO^l + 8x y yc^^) 

Proceeding to take derivatives with respect to Rj^ we have 

o 

^MO 1 ^ 1 (^l) 

d 


(5.57) 


d 

dR» 


where 
f 


O 

" XivyinM, (4x^y + — ) 

'"MO 1 \ o„ /?5 — !?/ 


2xy^^' 


The next derivative is 


.^ .^(,)(^ >2 = a {X]^o(^lb'*’®* yy^D ^^l)"^3 (^l)} 

dR^ o 

O 

^Mo (yib ®*^y ^ ®“i) 


(5.58) 


(5.59) 


(5.60) 


where 
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^3(^l) * ■ ^0 * (^Ib (^0 ■*■ ^^o)] 

o 

^ ^l^MO (^Ib ■*“ y VCjj k) (5. 61) 

o 

Two more derivatives yield the following expression 

'•'■u'2 * “ ■ (% ’■““l) • 

dRL o 

- -2'1 («l) - f'l' («l) 

o 

■" ^Mo (nb W'" ^'""i) <5* ®2) 

o 

Finally, this can be rewritten as 

+ af "* + af" 

+ u" [2a - (2Z yr )^1 + a^Uj = 1 ? ? (5. 63) 

^MO (^Ib ®*^y 

o 

The homogeneous part of this equation is seen to be identical to the 
equation satisfied by the optimal Uq. The forcing terms are seen to be 
functions of the zero-order optimal state and control. 

Having derived equations for the optimal Uq and u^ in the left side 
boundary layer it is easy to do the same for the right side boundary layer. 

Since the development is a direct parallel of that of the left side boundary 
layer, the details will be omitted. The results will be presented, 
preceeded by some comments regarding differences that appear in the 
right side results relative to those of the left side. 

First, because of the assumption in the equation for M that M ^ 0 

and Q 1 in the boundary layer, we define M and Q as follows 
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M 


- m 



~ ni. 

M' » -i M' 


“j5(h - h, + Imx) 

Q = e ^ 



^ 5 . 64 ) 

( 5 . 65 ) 

( 5 . 66 ) 


That is, mass Is referred to its final value and altitude to the value 
that would obtain on a Breguet cruise when mass reaches its final value. 
As a result of the definition of M, we have the costate relationship 


f 

X— = — ^ \ 

m^ 


( 5 . 67 ) 


Second, because of the stretching transformation 


Rj - R 


( 5 . 68 ) 


and the resulting derivative relation 


d _ Id 

HIT ■ ■ c 


( 5 . 69 ) 


all three state equations will have the signs of their derivative terms 
reversed. This sign reversal will, of course, appear in the Hamiltonian 
in the inner product of the costate vector with the differential equation 
of the state vector. As a result, all three costate equations will also 
have the signs of their derivative ten.as reversed. Third, since boundary 
conditions on M and Q are r*«w specified at R - R^, the constants of 
integration in the state and costate equations will be different from the 
corresponding values from the left side. 


Despite these differences, the equation for the optimal Uq in the 
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right boutidary layer is identical to that of the left. 


Uq + u"q (2a - [2Z JK]^) + a^UQ = 0 


(5. 70) 


The equation for the optimal is also unchanged, 

- af + af " 

I V 2 2 12^ 

u^ + u"j (2a - [2Zva]'^) + a'^Uj = i = i- 


(5. 71) 


3 


o 


') 


The apparent difference in the sign of the forcing term af *2 
fact that odd-power derivatives with respect to have the opposite sign 
from the corresponding derivatives with respect to . Expressing all 
derivatives with respect to R, no sign differences cccui’. Finally, it 
should be noted that the coefficients and the forcing terms in the above 
equations are defined exactly as they were in the left side boundary layer. 


5. 4 Matching the Optimal Initial Transient to Cruise 

We proceed now to investigate the solution of the differential 
equation for the optimal Uq and to determine the conditions under which 
solutions in the boundary layers can be matched to the cruise solution. 
The general solution for Uq on the oscillatory branch ^Fig. 5.1) is 

Uq = (Uqj cos u: Rj^ + Uq 2 sin a; Rj^) e 

+ (Uq2 cos tc Rj^ + Uq^ sin tc Rj^) e ' (5, 72) 


where and |J, the undamped natural frequency and the damping ratio, 
are given in Eq. (5.49), (5.52), (5,53). The frequency of oscillation is 



(5, 73) 
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If matching is to be possible, then the outer expansion of the 
boundary layer solution must be finite for small values of f . This means 
that the solution can not have positive exponential terms. Therefore, 

Uq 2 and Uq^ must be zero. Now we proceed to take the outer expansion of 
Uq from the boundary layer. First we transform the independent variable 


R 


L 



(5. 74) 


and then we take the limit as c -♦ 0. It is seen that 


Urn e 
<-♦0 


-Cu; R/€ 


= 0 


(5. 75) 


and hence the outer expansion of Uq from the left boundary layer is 
simply zero. This will be written in notation similar to that of O'Malley 



(5. 76) 


The use of the superscript i denotes the inner solution associated 
with the initial boundary layer. The superscript f will signify the inner 
solution associated with the final boundary layer. The outer solution 
represents cruise and is identified by the superscript o. 


At this point the inner expansion into the left boundary layer of all 
of the state and costate variables from cruise will be determined to first 
order in c. Obviously, no transformation is necessary to determine inner 
expansions of variables that are constant in cruise. These include Uq, u^, 

Vq' ^ 1 ' ^0’ ^MO' ^vO ^HO’ The others require transforming 

the independent variable and applying the limiting processes described in 
Chapter II. We have 
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(5.77) 


[Afl“2(.RL)] ■ (ri<R, 


+ €• 


2xyj (1 - x"^) 


- c 




from which 


2xJT^ + 6x^-3x^)- 

o o 


f . Ao2,i2 , , p . 2xy^b (1 - X ) , 

^ ^ yCp K (I + 6x^ - 3x^) ^ 


(AHq'^I' = 0 


A o,i 


2xy,. (1 - X ) 

(AH"1‘ = y,. R, 2 3- 

^ ^ yr^a + 6x^ - 3x^) 


Also 


[M^^^icRj^)] = (l - e 
from which 


^ ^2_“'yib^L^ / ^Ib 


+ € e 


^2x 

o 


(Mo°]‘ = 0 


^ ■ ^Ib^L 


Finally, 




which yields 


1^°'* ■ ><1. 


These values are summarized in Table 5.1. 


Now Uq from the left side boundary layer has been shown 


(5.78) 

(5. 79) 

(5.80) 

(5.81) 

(5.82) 
83) 

(5.84) 

(5.85) 
(5. 86) 

to have 
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Table 5.1 INNER EXPANSIONS OF CRUISE VARIABLES INTO THE 
LEFT SIDE BOUNDARY LAYER 


Cruise 

Variable 

Zero Order 
Expansion 

First Order 
Expansion 

y 

0 

■^Ib 

AH 

0 

2xy,. (1- x^) 

XuRt = n t-+M,^(0) 

0 

M 

0 

>lb^L ^lb^°> 

u 

0 

0 

S 

0 

_ 1 
~ a 


0 

1 

2 X VCp K 
^o 


1 

M., ^0) 
R 

^Ib 

1- ’lb 
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an outer expansion of zero, thus matching the inner expansion of Uq from 

cruise. We can write the boundary layer solution as 

-CUJ R, 

. n L 


Uq = (Uq^ cos cc Rj^ + Uq 2 sin w Rj^) e 


(5.87) 


In this forn», integrals of Uq will have the same form as Uq, that is 


^0 " ^^01 ^ ^02 ^ ® ^03 


(5.88) 


Obviously, in order to match Vq = 0 from cruise, is required to be 
zero. Furthermore, 

AHq = (AHq^ cos w R^ + AHq 2 sin o) R^) e " L ^ ^5 ggj 

A 

and AHq 2 = 0 to match cruise conditions. 

From inspection of the costate equations in the boundary layer it is 
obvious that Xj^q, X^q and Xjjq are all constrants. Matching to cruise 
conditions is simply a matter of equating these constant values to the 
corresponding cruise conditions shown in Table 5.1. The values required 
of XyQ and Xjjq agree with those values needed to minimize the zero-order 
Hamiltonian in the boundary layer. The first order Hamiltonian must also 
be zero everywhere along an optimal trajectory. If we consider its outer 
expansion for small values of c we can take advantage ^f the fact that the 

A 

outer expansions of Vq, Uq and AHg are all zero. The first order 
Hamiltonian then becomes 


^•^1^1° = ■ ^ ^MO ^Ib " ° 


(5.90) 


which iriplies that 
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(5.91) 




^MO ^ y. 


1_ 

'lb 


agreeing with our matching value. 


Three first order variables in the left side boundary layer are 
completely determined by the zero order state and optimal control. 

These are M^, and Xj^|^. Consider first the costates. 

o 

Every term of the integrand is of the form of Uq, a damped sinusoid, and 

hence integrates to the same form. Thus the outer expansion of Xjjj is 

simply the constant of integration. For matching with cruise, this constant 

must be , from Table 5.1. 

2xyc' K 
^o 

The solution of the canonical equation for Xj^^ contains one secular 
term, and all other terms are damped sinusoids. 

° ^MO S * “o y 

^ (’■O ^ ^0)} «L '5.93) 

The outer expansion of (X|^''l is given by 

(Xmi 1 = Xjyjo ^ib 

The coefficient of Rj^ is equal to one, because of the value previously 
assigned to Xj^q, and matching occurs if the constant of integration is 

(Mj^,(0)/yjb>. 
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Now consider the equation for Mj 

“i ■ Kb (? *“o *T~==) * yd‘ (tt * 


V / * V 2 




+ (4»Vyr^+4.)(u„ + AH,) +4x^v„(u„+AH„)}dRL (5.95) 

o 

Again, all terms are damped sinusoids except for one secular term. The 
outer expansion of is 

= Vjb Rl ®i 


where c., the constant of integration, has a value such that the initial 
value of Mj is zero. Thus the specified initial value of M is wholly 
satisfied by the initial value of Mg in the boundary layer. Now if M from 
the boundary layer Is to match with M from cruise (Table 5.1), then c^ 
must equal Mj^(O). 


If c^ happens to be zero, then the initial conditions that determine c^ 
will describe a locus of points from which transition to cruise can be 
accomplished at the same schedule of range and mass that would be 
experienced along a pure Breguet cruise. This locus will, of course, 

A 

include the origin in Vg - AHg space. In general, Cj^ will not be zero and 
then Mjj^(O) represents the initial mass of the Breguet cruise that matches 
the transition. A negative Mjj^(O) represents the fuel penalty incurred in 
climbing to cruise and a positive represents a fuel saving as, for 

example, in transition to cruise from a higher altitude than cruising altitude. 


The first order control, Uj, will have the same form as the zero order 

control, namely a damped sinusoid, but it also has forcing terms. The 
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f 


.J 


outer expansions of the forcing terms must be determined before Uj can 
be matched. 


IV 

The first of these forcing terms Is f^ where 

'i ' ■ ^Mi [nb ’’o ("o 


+ ^Mo“i (2yib * 

o 

Notice first that the coefficient of is a constant. Also, from Table 5.1, 
the inner expansion of is a constant plus a secular term. Only two 
derivatives of will remove the effect of its secular term, and hence 

IV 

the contribution of to the outer expansion of f^ is zero. The same is 
true of the term The remaining terms involve a sum of damped 

sinusoids (within the brackets) multiplied by which is itself a damped 

sinusoid plus a secular term. Every term in this product will be multiplied 
by a decaying exponential and will have an outer expansion of zero. The 
other terms, fg and fg, are of the same form as f^ and since each of them 
is differentiated at least twice in the forcing function, they too will 
contribute outer expansions of zero to the forcing functions. Thus the 
outer expansion of Uj^ will not be affected by the forcing functions. As was 
the case with Uq, the outer expansion of u^ will be a constant and the value 
of that constant will be zero, to match the value in cruise. 


Integrating u^ gives us whose outer expansion will be a constant. 
The value of this constant is seen from Table 5.1 to be the first-order 

A 

flight path angle from cruise. Finally, integrating gives us which 
will have a secular term from the integral of 
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(5.98) 


I 


IAH,r=y,^RL»c 

The constant is selected from Table 5.1, and matching of state, costate 
and control between the left side boundary layer and cruise is completed. 

We now briefly consider the optimality condition that = 0 over 
the entire optimal trajectory. Considering those variables whose outer 
expansions are non-zero we have 

(.X^2 1° * ^M1 ^Ib ^MO ^Ib (■ ^1 2 k ^ ^ 

^o 

Substituting outer expansions we have 


"HI 


1 

2k 

o 


(5.101) 


which is consistent with our previously determined matching conditions 
on 


5 . 5 M atching the Optimal Final Transient to Cruise 

Matching conditions at the right side boundary are similarly 
established. The development is a direct parallel of that of the left side 
boundary b*. t with different values required for matching. Table 5.2 shows 
the inner expansions of state, control and costate variables from Breguet 
cruise when extended into the right side boundary layer. These values 
were developed by stretching the independent variable in the vicinity of H 
equals Rp for the Breguet solutions by means of the transformation 


( 
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Table 5. 2 INNER EXPANSIONS OF CRUISE VARIABLES INTO 
THE RIGHT SIDE BOUNDARY LAYER 


Cruise Zero Order 

Variable Expansion 


First Order 
Expansion 


y 




'lb 


ViK ^ f n 

y (R -R )+—i2 — £2. +M,.(0) 


2xy^ (1-x^) 

VCj 5 K (l + 6x^ -3x^ 


M 


u 


1-e 


-y,uR 


lb fo 


0 


' ^‘"""’[nb(Rn-«R) ^ “ib«» 

^ ^Ib ^fO 1 

2xyU^-' 

o 

0 


X 


y 


0 


\ 

a 


X 


H 


X 


M 


0 


1 .’'lb ^fO 

’ib 


2x^/CJ^ K 


.’'ib^for. 


^ ’'lb ^fO 

LRj^.Rr+-^ + 


'^Ib 


2x,/e^ 

o 
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R = - cRp 


(5.102) 


and expanding the dependent variables for small c. Also, since is 
specified, we have at Rj. 

Myj. = (5.103) 

Mjj. = 0 (5.104) 

and the unspecified value of R^. can be expressed as an expansion 

Rf = Rof + cRff (5.105) 

thus accounting for the fact that since the first order solution does not 
change the final value of mass, it must change the final value of range. 


In the right side boundary layer the optimal control on the oscillatory 
branch has the form 




(Uqj cos u; Rp + Uq2 sin tc Rp) 


+ e 




(Uq 3 cos tc Rp + Uq^ sin ui Rp) 


(5.106) 


The coefficients in the divergent term are taken as zero so that matching 
will be possible. The exponential term is transcendentally small for 
small c 


0 


(5.107) 


It is seen that the outer expansion of u^ from the rignt side boundary layer 
is zero, which matches the constant value of zero from cruise. 
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It follows directly that in the boundary layer is a damped 
sinusoid and its constant of integration is zero. Furthermore, in 
the boundary layer is also a damped sinusoid and its constant of integration 
is also zero. That is, from the definition of AHq, namely 


AHq ' AHq - l»X - Rq^ 


(5.108) 


we will have at AHq equals zero 


AHq - lux + Rqj 


(5.109) 


as shown in Table 5.2. 


The form of the equation for M is unchanged from the left side 
boundary. It is again a secular term plus a constant of integration. All 
secular terms in expansions from the right side boundary layer are 
oppos te in sign from their counterparts in the left side boundary layer 

because of the previously noted sign difference in the differentials dRj^ 

» 

and dRp. Continuing with M we have 


= - y^^ Rp + c^ 


(5.110) 


" Or) nb^R * ^^f 


(5.111) 


In order to achieve a specified value of we must have the 
constant value of Mq equal to M^. and the constant, c^, must be chosen 
such that the value of is zero when Rp equals R^. We have 


(m/2) . m, + c [ - g) 


’'lb®R * 't] 


(5.112) 


This must match the expansion from cruise which is, from Table 5.2, 


V 





2xvq7x 


(5.113) 


Now we have the following condltioni’ ‘or matching 


(5.114) 


'^Ib^Of 


(5.115) 


Cj. = e 


■VlbRof [ * M,(0) . Vj,R„] 


•2xyr^ 

o 


(5.116) 


The first condition defines Rq^ and is in fact a zero-order statement 
of the Breguet range equation. The second condition is equivalent to the 
first since by definition of Mj. 

■ViK^nf 

1 - e = M, = 1 £ 5.117) 

f m, 

The third condition relates c, and R^j.. Since we know that equals 
zero, the equation for Mj evaluated at R|. serves to define c in terms of 
the zero order terminal state (AHqj., yQf). The third condition therefore 
defines Rj^. in terms of Cj.. The derivatiion of the required value of c^ is 
deferred until after matching has been established for the remaining 
variables. 

The final value of mass on the matching Breguet cruise is 


M,° = l 


[MjblOl 


^ ^Ib ^Of . . oil 

2xyrrTT 


(5.118) 




(5.119) 


V 



which does not in general equal the specified final value of mass. The 
difference is the first order term above and it is attributable to three 
sources. First, the initial value of mass on the matching Breguet cruise 
may have a non*zero value, Second is the fuel penalty paid to 

account for the climb at constant flight path angle that is characteristic of 

V 

the Breguet cruise, ■ . - ,. i: . Third is the penalty p. '.d in order to 


achieve a terminal state that may not be on the Breguet cruise, VnjRjf 
This is discussed more fully at the end of this chapter, where a more 
general range equation is developed. 


As before, the costates and Xjjq are constants whose values are 
seen to be both zero. The value of Xj^^jq, another constant, is determined 


from the requirement that.X'j = 0 

] = - 1 + XJJq VjJj = 0 


(5.120) 


^MO = 


(5.121) 


But since 


^MO ' mT ^MO 


^5.122) 


we have 


, .”i I - 1 


MO m^ 


(5.123) 


The solution of the equation for X^| again involves a damped sinusoid 

and a constant of integration. The constant is chosen as - and 

2xytj3 K 

o 

matching is achiev«.d for Xui. The equation for X«, again yields a secular 


V 





term plus a constant of integration 


fiO 


] = - X^o ^Ib ^ 


r\ fiO 


^Ib^f 




(5.124) 

'5.125) 


The value of the constant is, from Table 5.2, 
/^'ib^Of r „ ^ ^ ■>'lb ^ 

r~ 

D 


[ R„ . V! t 

*" ^Ib 2 x7C7,"R-' 


The form of the forcing functions in the right side boundary layer 
is the same as it was in the left side. Again, they contribute nothing to 
the outer expansion of u^. The constant of integretion associated with u^ 
must be zero. Now has only a constant as its outer expansion and that 
constant must be The non-zero value of introduces a secular 

term into the expression for We have 


«- p 


■ - nb>'R " 


and the value of c is 

R«f 2xy,. (1-x^) 

c » - 1 ^ 9i — 2 2- (0)+ y,. R,. 

2x/(!:”K K (l + 6x^-3x^) 


(5.126) 


(5.127) 


The optimality conditions on the Hamiltonian are consistent with 
the matching values of these variables. For = 0 we have 


'■yO “ ^HO “ ° 


(5.128) 


From.Ji'^j = 0 we have already established that 

>"lb^0f 


X .(-L'l 

^MO Wj 


(5.129) 


^Ib' 
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For .^2 = 0 we have 


(Xj.,1 1° = - - X^Q [- Mj + —j 

o 

,X f|0--rp -P + + ’'lb ^Of 1 

''hi' LHu 


(5.130) 


rihO^ir-RBl-MlKHn- 

’'ib '-2X yr]7Fr H R lb 

O 

1 


y 

'lb ”0f 


2x^Cjj K 


;] (5.131) 


'•m 


2x K 

o 


(5.1.32) 


which is the vMue required for matching, from Table 5.2. 


5. 6 Composite Matched Asymptotic H^xpansions and Costs in Transitions 

Having established the conditions for matching cruise to both 
boundary layers it is now possible to express the optimal values of state, 
costate and control variables in matched asymptotic expansions. These 
will be be uniformly valid over the optimal trajectory between the point 
at which the trajectory leaves the constraint of maximum Cp and the 
point at which it meets the constraint of minimum SFC. 


For this problem a matched asymptotic expansion of a variable will 
consist of the sum of the solutions for that variable in cruise as well as 
in the two boundary layers, minus the inner expansions of this variable 
as it passes from cruise into the boundary layers. Using the notation of 
O'Malley we have, for mass 


..c2 ..i2 , ..o2 . ,.f2 

M = M + M + M 




(5.133) 
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In order to proceed with the expansion it is necessary first to solve 
the zero-c.'der state differential equations (first order for M) in the 
boundary layers. Consider the expression for the optimal value of Uq 
in the left side boundary layer as a damped sinusoid 

Rr 

. n L 


Uq = (Uqj cos oj Rj_^ + Uq 2 sin co Rj_^) e 


(5.134) 


Integrating twice and using the previously identified constants of integration 




which can also be written as 


>0 


"C‘*^n^L f '' 

= e Ivqj cos uj Rl + Yq2 sin a; R^/ 


(5.135) 


(5.136) 


and 


A e 


or 




n 


7 — {■ («“ nV * “->'02)'=‘>® “ "l * t“n’'02) “”l} 

(5.137) 


AHQ = e 


-Cw R 


cos uj Rj^ + AHq 2 sin U3 R^} (5.138) 


The initial conditions on and AHq are obviously and 

A 

The remaining constants, Uq^, Uq 2 » ^02 ^°nld be expressed in 

A 

terms of and AH^^. However, a more useful relationship exists 

A 

among Uq, y^ and AH^. Consider the following definitions 


>01 ■ • «“ n “01 ■ ““02 


(5.139) 
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>02 * ■""“0r«“n"02-* “Ol ' (s)(’'02 ’■<“n"o2) 

(5.140) 

^"Ol “ • (t“n>'01 * “’' 02 ) 

(5.141) 

"n 


^**02 ' ( 7 ^) (“’’01 ■ **“n’'02) 

(5.142) 


Substituting for and then Uq^ from (5.139) and (5.140) into (5.142) 
yeilds 




^*^02 '^02 * “ ^ ( i ) ^02 

CO 

n 


(5.143) 


Substituting (5.139) and (5.140) into (5.141) yeilds 

^^01 “oi ' ■ ^ (“t) >01 

(tl 


(5.144) 


and from these two results we have the composite result 


AHj + tl, . - 2 (i) 


(5.145) 


Now simply by regrouping terms, the M equation can be written as 




’'0 


2x7U^ 




l(“o * 


2x^Cjy k 


(5.146) 


Making use of our result for Uq + AHq this becomes 


M. 




>■0 


o 


:) 
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The first term is easily integrated by relating Uq and yg to the 
differential equations for Vq and AHq respectively to give 


A, y.' AH' 




2x^Cq k' 


2xv^Cj^ K 


To integrate the second term, express as a damped sinusoid 

-CWnRi N 

^0 ' ® 1^01 ^02 ^ ^ l ) 

> 0 ^ = e ^ (yjjj2 cos^ uiB^+ sir? oj 


2^01 ^02 ^ “ ^l) 

-2Cu) R. 

S>'o^‘‘^L • i (“V • “’’i 


(5.150) 


02 


n 


^C^n^Ol ^02) “ ‘^o® 2 wRl (Cw„yoi^ V02^-^2ajyQj y^^)] 


(5.151) 


Summarizing 

Vq ^Ib 




2xyU^ 


n 


O 

+ sin 2a, Rl (a,y„,“ - - 2Ca,„y(„ y„j) 

- cos 2a, Hj_ (ca,„yji^ + Co,„yo2“ * ^a,^^ y^^)] + Cj <5. 152) 

The coefficient of the exponential is simplified by substituting for 

2 

U from its definition (5.49), (5.53) and then combining like terms. Then 
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after substituting the appropriate functions of for and AHq, and 

A 

substituting for ^ terms of and AHq^ we have 

^Ib ”^^n^Lr 1 / 2 ^ \ 1 

1 = ' [’’0l‘=““”L-3("n “”*lJ 


M, 




2xyU^~K 


[AH„,cosa;Rj_+i (y^j + C<c„AH„,) sin a- R J 




4u,' 




+ u: sin 2<s Ej_ (y„j^ - AH„j) 


2( jc ^ 


* C»^„ cos 2s, Rl (y„,^ + -j2 AH„,y„j + + Cj (5. 153) 

Since the initial value of Mj must be zero in the boundary layer 
solution, we have for the constant of integration 


‘=i“’'ib 


'Oi 


A 

ilH 


Oi 


“ ^o ^o (5.154) 


A similar expression in terms of Rp, and applies in the 
right side boundary layer. The signs of the expressions for y^j, yQg. 
and AHgj and AHq 2 opposite from those of the left side boundary 

layer. This leads to a change in sign in the expression relating the 
optimal control and the state 


"o*^o'*2(i)y, 


(5.155) 


The sign reversals are due to the reversal in sign of the derivatives 
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of the state variables in the right side boundary layer. We have for 

“i ' ■ nb"R “"r ■"^(v “ot ■ ?“n’'ot) " '*rJ 




■2 C“„«r , 
e L 


[aHojCos u^Rp sin uiRj^] - £— j 




^ V ui n ui ' n ui ui/ n \ ui 

^O 

2 oc 

- y Cn;„<=os2n.R„ {yj - -f V) j 

From this we can write 

Mj'e '{-yibRp+^e ' " ”L'>'0t=°="®R*3(‘Hi^*«0ft“n%)“‘"‘"'^R] 


■ ’'Of) =*" “ "r] 

■^Cu.' Rr 

■ (’'lb T (’'O 

o 

2 cc 

+ OJ sin 2«R JI (yj - t Ci«„ cos 2a.R„ (y^j^ ' “T *"of ’'of 




(5.157) 


The value of the constant is such that is zero at the final value of 
R (when Rp equals zero). The specified final value of M will be satisfied 
exactly by the final value of Mq. Accordingly we have 
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AH 


Of 


< 2xy^ 


) 


* ?l53;(’'lb * (’'Of^ * "n^^®0t^)] 

o 

when the final state is fully specified. 


(5.158) 


Transforming both boundary layers to functions of R we can write 
out the constituent parts of the composite solution for M 
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(5.160) 
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( 5 . 163 ) 


Finally, the composite solution for M to two terms is 
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^ ^02 ^ . |-j^o2Jl2 . ^^o2J2 

At this point we can note that similar composite solutions can be 
expressed for u, y, AH, and the costates, and the matching constants 
have already been evaluated. In all these cases, however, the results 
will depend on u^. While u^ is, in principle, evaluated quite directly, in 
fact the high order derivatives involved in the forcing functions make an 
analytic evaluation quite laborious. It is shown in Chapter VI that a very 
good representation of the optimal trajectories comes from considering 
a zero order boundary layer solution to match its corresponding cruise 
solution through first order in c, and the labor involved in evaluating the 
first order corrections to the transients is not justified. 

12 

We have seen that the constants of integration associatedwith M 
f2 

and M are an indication of the difference between the fuel consumed in 
transition to and from cruise and the fuel that would have been used in a 
pure Breguet cruise over the same interval of range. The equations 
defining them have the same region of applicability in AHq * Vq space as 
do the state equations, and matching between cruise and transition to and 
from cruise is possible everywhere within that .region. It is possible 
therefore, to assign a cost number for any point (AHq, y^) in the region 
and to develop contours of constant cost. Each contour will define a 
locus of initial conditions from which the same cost is incurred in 
traversing a matching transition to cruise. The same could be done for 
final conditions on transitions from cruise. The cost is easily evaluated 
relative to the cost incurred in a pure Breguet cruise of the same range 
using equations (5.154) and (5.158). 
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The contours are not optimal trajectories. Trajectories, as we 
have seen, are damped sinusoids as a function of R. In AHq - Vq space 
these become spirals. In traversing a spiral one passes through a region 
of higher as well as lower initial cost. This is because the spiral 
trajectory may include a region of n jative AH wherein the airplane is 
climbing and requiring more thrust and hence a higher fuel consumption 
than for cruise, as well as a region of positive AH wherein the airplane 
could fly down to cruise at reduced thrust. 

Since any point in AHq - space within the region of applicability 
of the equations can lie on an optimal transition trajectory, it is possible 
to assign to each point the throttle setting or thrust coefficent that the 
optimal trajectory would requi.e as it passed through that point. Obviously, 
all optimal trajectories will include the origin, at which point the throttle 
setting equals that required for cruise. Moving away in one direction all 
throttle settings will ultimately reach maximum. In another direction all 
settings will reach minimum. These and other loci of constant thrust 
coefficient are discussed in Chapter VI in connection with some numerical 
examples. 


5. 7 Non-Oscillatory Optimal Control 

Thus far we have considered only the oscillatory form of the 
optimal control. It is, of course, possible that the roots of the 
characteristic equation of the differential equation for the optimal control 
will all be real. We then have 


+r,R ■•■‘'2^ 

^0 ' “01 ® ^ ^^02 ® “03 ® ”04 ® 


(5.165) 
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The roots are symmetrically located about the imaginary axis so that 
we have 


*'1 “ ■ **3 


^2 “ ■ *-4 


(5.167) 


The development of the boundary layer solutions and the matching 
of these solutions to cruise proceeds In almost Identical fashion to the 
oscillatory case. If there Is to be matching, the coefficients of the two 
divergent exponentials (say Uqj and Uq^) must be zero. Then, by directly 

A 

Integrating the zero order differential equations for and we have 


•“Ol /■'l” ^ “02 
— e + e 


r“oi 1' 0 


(5.168) 


AH - + a r^Ol .‘*■*'1*^ . “02 

AHq - + a L — 2 ® — 2 ® 


(5.169) 


We would like to express the optimal control, u^, as a linear 

A 

combination of the two elements of the zero order state, and AH^. To 

r R r^R ^ 

do so we solve the above two equations for Ug^e ^ and u^g ® in terms 

A 

of Yq and AHq and then express their sum, Uq, as 


r, +r. 


“0 “ ■ (^) * (\^) 


(5.17U) 


Th' s Is the optimal control for initial transition. For final 
transition (away from cruise) the control law Is slightly different due to 
the change In sign associated with the state equations In the right side 
boundary layer. The control law for transition from cruise Is 


r, +r. 


“0 • - (^) - (V^) r, 


(5.171) 


r 


In order to derive an analytic expression for the cost in fuel 
associated with transition we consider the state equations in the form 


rjR 


^0 "^ 01 ® >02 ® 


^2^ 


A A ^1^ A ^2^ 


0.172) 

(5.173) 


Initial conditions are defined as 


''Oi 


A 

AH, 


>01 >02 
01 “ ^^01 ^^02 


(5.174) 

(5.175) 


Now from a comparison of equations (5.168) and (5.172) we can 
write for equation (5.174) 


>0i “ ^ %1 '‘02 


(5.176) 


Similarly, equations (5.169) and (5.173) allow us to write for equation 
(5.175) 


^01 * “T '^01 ^*02 

r, r« 


(5.177) 

Ufii 

‘1 ‘2 

The right side terms of equation (5.176) will change sign for descent. 


Now equations (5.176) and (5.177) can be solved to express u^^ and Uq 2 

A 

in terms of the Initial conditions, and AHq^. 


.. , X /'*'A (’'or*'2^*^oi) 

,, . X ('“l^”oi">Ol) 


(5.178) 


(5.179) 
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For descent, equation (5.179) Is unchanged but equation (5.178) becomes 

(5.130) 




From a comparison of equation (5.172) with equation (5.168) we have 


^01'* 


/ a 


^1 

(5.181) 

“02 

(5.182) 


with the negative signs applying in final transition (descent), and from 
equation (5.169) and (5.173) we have 


^^01 * (A) “oi 
*^1 

^^02 ' ( A )“02 
*^2 


(5.183) 

(5.184) 


The fuel cost associated with the initial transition is expressed by 
the equation 


D. 




4x^Cp K 2X4/Cp K 


L 

(5.185) 


For descent 
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axTCp; 


4x^Cj. K 2x,/C^ 
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(“o-^^oB^r] 


(5.186) 


To evaluate these equations we need the following expressions 

(5.187) 


2 .^**2^ 


n n2r«R (r.4’rn)R 

96 


p > p > ( 0 > 


(5.188) 
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Now define the following constants 
0 • ^ (“Ol * '^"0l)("02* ^®02)(nb" Sx’yvrs^- Xoi’'02(nb* 

♦ 4 x ^ [(“ 01 * Ai )’' 02 *(," 02 *‘*” 02 ) v ]} 
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* 4 xV ( uo ,* AHo ,)]} 

2>{(uo2*Afto2) ("r *^”^5'-/^^)*>'02[("r V 
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•► 4 xV ( U(,j + AHjj )]} 


(5.189) 

(5.190) 

(5.191) 

2 )R 

(5.192) 

I 

(5.193) 

(5.194) 

(5.195) 


97 


The ascent cost can now be written as 


M, = y,uRT 


^ib 


^ L ^ ““0 ■ ^2 

a^ 2 tjRj^ ag 2r2RL 
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(5.196) 


where 


3> y 

a ^Oi 


A A A 

^Ib _ ^0 ^1 ^2 
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(5.197) 


The corresponding cost for descent transition is 


Mj = e 


■^Ib^f 


y R +!lb ^b 

.'Ib^R '0 


2^C:^ Kx 
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“O ,*''l'''‘'2*’*R . “l , “2 ,^*'2®R 


] + C^ 


(5.198) 


Cf = e 


^ ^Ib g +_l0_ + ^+^l (5 199) 




5. 8 Solution Without a Cruise Section 


If the specified final value of mass is sufficiently large the solution 
will not achieve cruise altitude. The end conditions on altitude and flight 
path angle may even lie within the initial boundary layer. For slightly 
smaller specified values of final mass there will be two boundary layers 
but they will coalesce, and no cruise section will exist. In the absence 
of a cruise section matching is not required. As a result it is necessary 
to use the most general form of the optimal control in evaluating a 
solution. Those constants associated with the positive exponential, which 
were required to be zero in c'der to achieve matching, may now be non- 


zero. 
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The optimal control, flight path angle, and altitude difference are 
now expressed as 

-«“n«L 


Ufl'e 


(Ujjj cos u; + Uq 2 sin a; Rj^) 


+C‘*;n^T 

+ e (Uq3 cos o; R^ + Uq^ sin a: R^) 


(5.200) 


>0 ' ® <^01 ^02 

+ e (Vqj cos (ti Rj^ + Vq^ sin u Rj^) 


(5.201) 


AHg = e ^ ^^^01 ^ ^02 ^ 

+^(0 R T A A 

+ e " ^ '*’ ^04 ^ 


(5.202) 


The twelve unknowns in these equations can all be identified in 
terms of the specified initial and final values of the state variables. 

A 

First, by differentiating the above expressions for AHg and yg with 
respect to Rj^ and equating the results to yg and aug respectively we 
have 


dAH, 


uttn^ r A A . 

- + C«„AHn,) sin u;Rj_] 


+ e 


01 ■ ’•‘-n“"'02 

+CuJnRL r A 
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n“"g4' — “"'l] 


+ (-tcAHgg + Cw„AHft,i) sin wR 


(5.203) 
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dyQ 

3 n ~ ■ ® L ^‘*^^02 ■ ^L 

- (wy^i + Co.'„yo2) sin a R J 

+ e l_(u'yo4+Ct«J„yo3)cosu.'RL 

+ (-ccyQ3 + C‘*;nyo4) sin u; R J (£ 

Now by comparing terms from these differentiated expressions 
with the original equations for and Uq we have 
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where the matrix A is defined as 
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We now have equations to express the eight elements of the vectors 

A 

and AHq in terms of the four elements of the vector Two more 
equations are available from the initial conditions 


>'i = ‘>^01 ^03 <^-210) 

AHi = A% + **Hq3 (5. 211) 


and two from the final conditions 
Vj- = e (y-Qj cos u; 

+Co;n^f 

® ^'^03 it Rf + 7Q4 ^ ^f^ 


(5.212) 





+ e 




(AHqj cos Or + AHq 2 sin u; R^) 

A A 

(AHq 3 cos u! R^ + ^04 “ ^f^ 


(5.213) 


These final conditions introduce another unknown, R^, and hence require 
another equation. The required equation is that of M^, the specified 
final value of M 


M, 


' W^L’'0-SoT) 


>0 


2VCj^ Kx' 
0 


* 4X^ J 


2VCp Kx' 


RL=Rf 


(5.214) 


This expression can be integrated directly. Some of the tedium is 
removed by referring to the manner in which this expression was 
integrated in the initial and final boundary layers for matching with 

A 

cruise. First, it is useful to compare the sum of Ug and AHg with yg 
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'"o ^ <y01®°®“=^L >'02 sint^’^L) 


n 


-C‘^n*^L 1 

■•■ ® <>'03 cosujRl + Vo4 sintcRL); 


(5.217) 


Now we can evaluate the combinations of products of the terms 

A 

Vq and Uq + AHq that appear in the Integrand in terms of the components 
of 2 q* We have 

2 ”^<*^n^L 2 

yp = e (y^j cosccR^ + 7 q 2 sinwR^) 

+2Ccc R, 2 

+ e " ^ (ygg coswRj^ + y^^ sinccRj^)'* 

+ 2(yoiCOSU)Rj^+ y^g sinwRj^)(yQ3 cosa'Rj^ + yQ^ sintc R^^) <5.218) 
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(Uq+AHq) = 4-^ |e (y^jCosu-RL ■*■ >02 sinaiPj^ 

+2Co;„Rl 2 

+ e (y^g cosicRj^ + sinu:Rj^) 


- 2(ypjCosccRj^+yQ2sina;Rj^)(yQ3COSu;Rj^+yQ^sina'h\ )j (5.219) 


A y -2CoJ R. 5 

= 2^ (y^j coswK^ + y^g sina:R^) 


- e (y^g cosu:Rj^ + sino-'R^) 


(5.220) 


Combining terms we can now express M as 
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-sX, + 4x^y 'f^)]} ‘"'l (5. 221) 


The following integrals can now be used, in the appropriate 
combinations, to express the integrated form of M 
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(5.222) 
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+ ojsin 2 ujH^ ± 2a)Rj^^ (5.223) 
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(5.224) 


4(0. 


r n 

1 — l=^t 
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R. sin 2a'R, 


± sin2(oR^ - cocos 2ujRj^J (5.225) 

^ I- cos^ dRj_ = -^ + ^ (5.226) 

^ ^ sin^ u)R|_ dRj^ 


Rj^ sin 2u;Rj^ 
= “2 ?w — 


^ sin ujRl cos coR^ dRL = " cos 2 coRl 


(5.227) 

(5.228) 


The two coefficients have previously been simplified as follows 

2^ (l + 4 -il l ± Sx^y ^ (l ±2x)fU^ X) = y +-~^± Sx^y X 


(5.229) 


Finally we have our equation for the specified value of in terms 
of the four components of and Rj^^ 
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3 ^^01^03 ■ '^02^04^ 3 ^^01^04 “^02^03^] (5-230) 


Since this equation is highly non-linear, an exact solution is 
impossible. An iterative solution can be obtained by guessing a value 
of Rj f, solving for values of the components of Vq, then evaluating Mj. 
from the above equation and repeating for a new gues.ied value of Rj j. 
until a satisfactory agreement with the specified value of Mj. is achieved. 


To evaluate the components of >-q we can easily simplify our set 
of eight equations in eight unknowns (Uq elements are not required). 


From our equations for AHqj and we can say 

>'02 ' " l^^Ol ^Ol) 

u: n 


(5.231) 
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Cw A 

** A TI 


‘^"04 ’03* -if ^03 


If we eliminate y^g and AH^g from the latter two equations by means 
of our Initial condition equations 


’03 ' ’01 ■ ’01 


^”03 = ^Ot - *"01 


and then substitute for Yq 2» ^^02 ^^04 ^ end-condition 

equations we have two equations in two unknowns 
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Now the solution proceeds as follows 

1. Guess a value of j. 

2. Solve equations (5. 237) and (5. 238) for and AHq^ 

3. Solve equations (5.235) and (5,236) for and AHqj 

4. Solve equations (5.231) and (5.233) for and y^^ 

5. Solve equation (5. 230) for and subtract the specified value 

of Mj.. Plot the difference against Rl|- If the difference is not 
sufficiently close to zero, return to step 1, guess a new value of Rj^j. 
and repeat until the difference at step 5 is sufficiently close to zero. 
Since the sign of dM/dRj^ must always be positive, the process of 
locating the zero crossing on the plot should not involve many trials. 


So far in this section we have considered solutions in which both the 
po'^.itively and negatively damped exponential terms make non-negligible 
contributions over the entire trajectory. If the specified final value of 
Mj. is steadily increased, the effects of the negative exponentials on the 
final state and of the positive exponentials on the initial state v/ill 
approach zero. There results 


^Oi ' ^01 


(5.239) 


^» 0 i = ^»01 


(5.240) 


.■^^‘‘'n^Lfi 


[y„3(cos«RLf 


t-Hof = e 




[AHp3(cosujRLf+-j2^8inu;RLf)+~sinu.RLj.] (5. 242) 
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The term with tne negative exponential factor now no longer appears in 
the equation for but otherwise that equation is unchanged. The 
solution proceeds in the same way as before, but the equations are some- 
what simpler. 


The equations of the state and optimal control now have the following 
form in the vicinity of the initial point 


Uq = e 


(Uqj cos + Ujj2 sin 

-Cw-Rt 

>0 ' e (yjjj cos wR^ + >02 

AHq = e " ^ (AHqj cos ccRl + aAq2 


(5.243) 

(5.244) 

(5.245) 


From the general expression for the sum of Uq and AHq we see that the 
control law in the vicinity of the initial point is 


"0 ' • ^"0 ■ s;; ’’o 


-i 


(5.246) 


which is identical to the control law in the left side boundary layer in the 
matched asymptotic problem. Similarly, in the vicinity of the final point 
we have 


+Cu.' R, 

^0 ' ® "*04 

+C« Rr 

>0 ' ® ^04 

A +C« Rf A A 

AHq = e ^ cos tcRj^ + AHq^ sin wR^) 


(5.247) 

(5.248) 

(5.249) 


Referring again to the general expression for the sum of Uq and uHq 
we find the control law to be 
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(5.250) 


which is identical to the control law in the right side boundary layer in 
the matched asymptotic case. 

The optimal control must pass from its initial form to its final form 
at some intermediate state. By inspection of the two forms it is obvious 
that the only requirement for that state is 

yQ = 0 (5.251) 

The altitude difference need not be zero since it has the same sign and 
coefficient in both forms of the control law. Thus in some short range 
problems the trajectory may not be required to reach cruise altitude. 

It is useful to compare this short range problem to the matched 
asymptotic problem. If in the short range problem the positive and 
negative exponential terms interact, it is as if the two boundary layers 
of the asymptotic problem were so close as to overlap. If the positive 
and negative exponential terms do not interact, it is as if the two boundary 
layers matched asymptotically to each other without an intermediate 
(cruise) section. 

Both forms of the short range problem are singular perturbation 
problems and require the transformation 

(5.252) 
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I 


to become regular. This transformation allowed investigation of an 
interval along the R-axis of width on the order of < and is consistent with 
our restriction to short ranges. Without the restriction to short range 
we must consider the cruise problem, which is of course singular. Thus 
by further extending the range, or final value of mass, we are led to the 
matched asymptotic problem with its initial and final boundary layers 
matching an intermediate solution representing cruise. 

5. 9 Breguet Range Including Corrections for Tran8ition^ 

In concluding this chapter we can evaluate an expression for final 
range. We have already modified the Breguet range equation to account 
for the gradual increase in altitude and the resulting increase in fuel 
consumption. We now can incorporate first order corrections to account 
for the possibilities that the initial and final values of and may 
not correspond to values on a Breguet cruise. 


From our terminal boundary layer matching condition we have 


1 - e 


"'^’ib^fO 



(5.253) 


That is, the specified value of determines the zero order value of the 
unknown final range 

Rfn » - -L i^(l - Mf) (5.254) 

Tib ^ 

The first order cor'n'ction to final range comes from the other matching 
condition 
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+VihI^ 
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The value of Cj is defined as a function of initial conditions on and AHq 
by (5.154) for C s 1 and by (5.197) for C > !• Likewise is defined as a 
function of final conditions on yp and AHq by (5. 158) for C ^ i 
(5.199) for > 1. We now have for the final range 


*^f ' ^fO 




(5.256) 
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The first term is seen to be identical to the Breguet range equation 
as it was developed in Chapter II; a correction factor appears, decreasing 
final range, to account for increasing altitude at constant fligiit path angle 
in cruise. The second term represents the change in range due to 
transition from an initiul state that is not on the Breguet cruise. It can 
be related to a non-zero value of mass on the matching Breguet cruise 
at equals zero. It can be expressed in terms of the initial values of 

A 

and AHq and is zero if they are also zero (i. e. , if they are on Breguet 
cruise). 


The thi .'d term represents the change in range due to transition to a 
terminal str.te that is not on a Breguet cruise. It can be related to the 
mass difference (at R equals R^) between specified final mass and final 
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mass on the matching Breguet cruise. It can be expressed in terms of 
the final values of and AHq and is zero if they lie on a Breguet cruise. 


Finally, by making use of Eq. (4.45), (4.81), (4. 83) and (5.253), 
the range equation (5.258) can be identified with the more familiar form 
of the Breguet range equation 
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o 

Both of the range correction terms can be positive, negative or zero, 
depending on whether the average thrust required for transition is less than, 
more than or equal to that required on a Breguet cruise over the same 
range. A positive value of c. will result in a reduced final range. This is 
to be expected, since a positive means that the Breguet cruise to which 
the initial transition matches has a fuel budget that is reduced at R equals 
zero by c c^ Similarly a negative value of c^. means that the matching 
Breguet cruise terminates at a value of mass that is less than the specified 
value, M^, and hence translates into a loss of range incurred in 
diverging from cruise to meet the specified final values of y^ and AHq. 

Both of these terms are zero if the initial and final state are on the 
Breguet cruise and in that case the range equation is identical to that 
derived from Breguet cruise in Chapter IV. 


112 





1^ r -A-v--?,* 




CHAPTER VI 

COMPUTATIONAL RESULTS 


6. 1 Introduction 

This Chapter gives the results of some computational studies of 
transitions and the costs associated with them. Several non- optimal 
control policies are described and comparisons are made between them 
and the optimal policy. The comparisons are made for three aircraft. 
The Boeing SST is used to represent aircraft that cruise at supersonic 
speed. The Boeing 707 represents aircraft that cruise at transonic 
speed. The McDonnell Douglas F-4 is used to study transonic cruise in 
an aircraft :hat is capable of supersonic flight. Finally, some comments 
are made about how the optimal policy might be implemented in a flight 
control system. 

6. 2 Alternate Control Policies 

It is useful to compare an optimal initial transition trajectory and its 
cost with a series of trajectories which use non-optimal controls. These 
transitions assume that the aircraft has already accelerated to its cruise 
speed and is attempting to reach cruise altitude and level off. The first 
non-optimal control assumes that the aircraft maintains maximum power 
setting until it reaches cruise altitude and then assumes its cruise power 
setting and levels off in zero time. This trajectory will be called the 
(Cplmax trajectory. In the second, the power setting is first set at its 
cruise value and the aircraft then eventually levels off at its cruise alti- 
tude. This trajectory 511 uc '‘•led the (C^) . trajectory. The 

X* CX U1S6 


U3 


third Is a constant rate of climb trajectory, identified as yj^Q^st* 
again the power setting and attitude are assumed to change instantaneous- 
ly to cruise values as cruise altitude is reached. The fourth trajectory 
is the optimal trajectory, 

All of the suboptimal controls represent a constraint on one variable 
in the equation for constant velocity. The exact equation may be written 
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( 6 . 1 ) 
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If Cp is a constant, whether or cruise other 

value, we must have for u 


u ■ - 1 ± SQRT 
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( 6 . 2 ) 


If y is a constant, then u must be zero and we have for Cp 
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If u is u^p^ for Initial transients 
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and we must have for C, 
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For the ^^p^const there is obviously a limit on the 

A 

values of and y such that the argument of the square root (6. 2) 
remains positive. This is easily identified in - y space. It 
corresponds to the condition 


u =* -1 


( 6 . 6 ) 


and from our definition of u this is equivalent to saying that lift is zero. 

A 

Obviously, at a particular value of AH there is a maximum value of y 
at which constant velocity flight can be maintained. That situation cor- 
responds to minimum drag since lift, and hence the induced component of 
drag, is zero. 


It is also true that in constant rate of climb trajectories there will 
be a maximum climb angle above which constant velocity flight cannot be 
maintained. This value is determined from considering the maximum 
value of Cp required for such flights. 

The maximum value cannot occur below cruise altitude, or it will be 
impossible to maintain the constant rate of climb, and the trajectory will 
become a (Cp^^ax trajectory. Accordingly, using and cruise 

altitude (AH ■ 0) we have 


^^F^max ■ 


o 


(6.7) 
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The optimal trajectory will be limited by a locus of points correspond- 
ing to a maximum value of Cp with u satisfied by the optimal control law. 

A similar locus will exist for the minimum value of Cp. 

Now using the zero order equations for altitude and flight path angle, 
(5.21) and (5.22), and the first order equation for mass, (5.23), it is 
possible to evaluate zero order trajectories and first order costs associ- 
ated with them for these four ontrol laws. The equations are integrated 
using a fourth order Runge-Kutta routine. Tlie atopping conditl: I'cr 

u^p^ and (C pJ^ruise attaining of steady state altitude. For the other 

two controls the stopping condition is the event of altitude exceeding its 
cruise value. 

These equations are the state equations for the linear- quadratic 
problem and small initial values are chosen so as not to violate the 
assumptions inherent in the equations. 

6.3 Boeing SST 

For our first cost comparison we choose an initial altitude of 

AHq. » -0.30 (6.8) 

For various values of initial flight path angle up to the respective maxima, 
trajectories have been calculated for the four controls and the resulting 
costs plotted in Fig. 6. 1. The airplane used in these calculations was 
the Boeing SST (Appendix A) in supersonic cruise. The variable plotted 
on the vertical scale of Fig. 6, 1 is the difference in the first order mass 
term between the indicated transition climb trajectory and a pure Breguet 
cruise of the same range, that is It is based on equation 

(5.23). 
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It Is seen that is a good approximation to minimum cost 

but its range of feasible initial flight path angles is severely restricted. 

It can be concluded that if ^ were identical to (C_) then 

F cruise F max 

these two trajectories would be identical to each other and their cost 
would be virtually the same as that for the u^p^ trajectory. 

In a narrow range of angles the y^Q^st also compare 

well with minimum cost. This comparison worsens as approaches its 
maximum feasible value and worsens rapidly as y^ approaches zero. 

For realistic attitudes associated with climb to cruise, the constant 
velocity transition must be either u^p^ or In general the recom- 

mended procedure for a pilot to follow in flying his transition to cruise is 
to accelerate and climb at until cruise speed is reached, then to 

climb at (C-r~) and constant speed until cruise altitude is reached, and 
then to level off at cruise altitude and speed in an unspecified manner [22, 
23, 24 ] . It is seen from Fig. 6. 1 that in this comparison the cost im- 
provement in terms of Mj ranges from 0. 38 to 0.20. To convert this 
number to a weight it is necessary to multiply by c and by the initial 
cruise weight. For the SST this converts to a weight of from 540 to 285 
pounds of fuel. 

Figure 6. 2 is a comparison of the zero order trajectories for the 
four controls In R - AH space. All start from an Initial flight path angle 
of -0.05 radians and an initial AHq of -0.30 or -6,240 feet. The opti- 
mal trajectory overshoots in AHq by 0.00127 which is equivalent to 26 

A 

feet. Figure 6. 3 shows the same trajectories in AHq - y space. This is 
essentially a phase plane, and optimal trajectories spiral into the origin. 







Figure 6,3. Altitude Transi^^nts in Space for \arious Controls 






Figure 6, 3 also shows a (C ) trajectory approaching from some 

r iTiaX 

much lower altitude than cruise. If pursued to its limit it would settle at 
the maximum cruise altitude of the aircraft. This trajectory acts as a 
separatrix for all other (C trajectories originating at other initial 

conditions. All other trajectories will fare smoothly into this 

separatrix and continue on to maximum altitude. None will cross it. This 
also applies to trajectories from higher altitudes than the maximum 
cruise altitude. They would fly down to the maximum cruise altitude re- 
maining on one side or the other of the extension of the trajectory from 
infinity (separatrix). The separatrix for ascent at is also 

shown. 

A 

For negative values of AHq the separatrix follows fairly closely a 
locus of zero lift at maximum Cp. This locus is also indicated in Fig. 

6. 3 It is evaluated by equating u to - 1. 0 in Eq (6. 2). Above this locus 
the flight path angle would be too steep to maintain constant velocity 
flight. 


Figure 6.3 also indicates a locus of points at which throttle setting 
is maximum if the optimal control is used. This locus comes from set- 
ting C„ to in Eq (6.5). Above this locus the throttle setting re- 

quired for a constant velocity range- optimal transition would be greater 
than the maximum throttle setting. If one follows the separatrix 

backwards to lower altitude, eventually the separatrix will be above the 
locus of maximum throttle for The point at which this intersection 

takes place is interesting because in ascending at from some 

large initial altitude difference this will be the point at which it is possible 



to begin using the optimal control. Figure 6. 4 is a sketch of that inter- 
section and an optimal trajectory from it. Approaching from some lower 
altitude at one follows the separatrlx (A - A') until point C. 

There one begins the optimal spiral Into 0, down- throttling all the way. 
Consider now point as an initial condition. The optimal spiral requires 
upthrottlln^ Initially and the throttle satiurates at B'. The extension of 
this optimal trajectory Is Indicated In a dashed line. From B' a 
arc fares Into the separatrlx and eventually comes out of saturation at, or 
very near, C, from which It follows an optimal spiral to 0. 

For the Boeing SST the separatrlx and the locus of at u^p^ 

essentially overlap In the vicinity of their Intersection. From studying a 
digital computer print-out of the trajectories in the vicinity of their Inter- 

A 

section, the point (-0. 745, -0. 300) In AHq - space was taken as the 
Intersection. Since this point would be well outside the linear- quadratic 
region, a comparison was made using the full state equations (5. 3, 5. 4, 
5.5) and the linear- quadratic optimal control (5. 137) which, for these 
equations, becomes a sub- optimal control. The comparative trajectories 
are shown in Fig. 6. 5. The optimal trajectory overshoots the Breguet 
cruise by a ah of about 0.03 (624 feet) and returns to meet the Breguet 
cruise at a value of of about 15. This corresponds to a range of about 
60 miles and would require about two minutes to complete. The difference 
In M between the two trajectories would be 0. 00076 which corresponds to 
a fuel weight savings of 487 pounds. This can be converted by the Breguet 
range equation to a range Improvement of 5. 44 n. ml. In Fig. 6. 6, u and 
Cp are plotted as functions of for the two transitions of Fig. 6. 5. 
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Figure 6. 4. An Optimal Trajectory with Throttle Limiting 
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ts for the Trajectory of Figure 6. 5 
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Figure 6. 7 shows in AHq - space a comparison of two complete 
climb- cruise- descent trajectories. Single arrowheads denote the u^^^ 

transitions and double arrowheads denote the (C_ ) (C„) . transitions. 

F max' F min 

The cruise segment (^l^axisfrom 0.0 to +0. 10) is conr.oon to both 
trajectories. Initial and final values of y^ are zero. Initial altitude is the 

A 

same as final altitude. The initial value of AHq is taken as -0.20 and the 
increase in AHq during cruise as t-o. 10. This means that the final value 
of AHq is -0. 30. These altitude values are kept small in order not to ex- 
ceed the assumptions inherent in the linear- quadratic problem. The in- 
crease in AHg during cruise can be related to a zero order final value of 
range through the constant cruise flight path angle and then to a final 
value of mass which must be the specified final value of mass. 

In this presentation one can see that the zero order ascent and des- 
cent and the first order cruise are the most significant parts of the trajec- 
tory. The first order corrections to ascent and descent would be of order 
c smaller and would not make an observable change in the figure. The 
zero order cruise, on the other hand, would be represented by the origin 
alone and would not fairly represent cruise. It is possible to speculate, 
therefore, that one could make a simpler approximation to the analytic 
representation of the solution by asymptotically matching the zero order 
boundary layer solutions to the first order cruise. 

We now proceed to evaluate some numeric results related to Fig. 

6. 7 in order to show the relationship to the Breguet range of the first 
order corrections to it due to initial and final transitions and to non- zero 
cruise flight path angle. The numeric values used for the SST flight 


parameters are shown in Appendix A. 
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The altitude difference re.. Ing from cruise is +0. 10. This 
implies a zero order range of 


0. 4649 


(6.9) 


and a final mass of 


Mqj = 1 - e a 0.09516 


The zero order range in nautical miles is 


range = 6^76 ^ ^ 715. 8 n. mi. 


The complete first order correction to range has been shown to be 


A + Jl- + .... 

vjb nb * 


The third term is the Breguet correction due to non- zero cruise flight 
path angle. Its value is 


= -0.9903 (6. 

grange = -.9903 < .^K ^ “ -3.4 n. ml. (6. 

Thus the range achieved on a pure Breguet cruise for which the final 
value of M is 0.09516 is 


range = 712. 4 n. mi. 


\ A: 




Since initial and final conditions are not on a Breguet cruise there 
will be increments in fuel or range associated with meeting the initial 
and final state. From digital solutions using the linear- quadratic problem 
we have 

u 


c. 

1 


opt 
0. 44124 


<^F>max 
0. 76852 


0. 53041 


0. 48846 


( 6 . 16 ) 


The first order correction to the problem due to ascent and 
descent transitions is 


Rjj » -2.0513 + 2.2312 = +0. 1799 


Ai*ange = +0. 6n. mi. 


(6.17) 
(6. 18) 


for a total of 7i3. 0 n. mi. For the (Cp)^j^j^ - problem the 

correction is 


ARjj = -3.5728 + 2.0548 = -1.5180 


(6. 19) 


Arange = -5.2 n. mi. 


( 6 . 20 ) 


for a total range of 707. 2 n, mi. and the saving of u . over (C_) 

opt X* 


max 




Looking at the components of the first order corrections for u . we 

opt 

see that the amount of range lost from a pure Breguet cruise because of 
optimal transition to cruise from a lower altitude is 


ar 


-c 


c. 

I 


( 6 . 21 ) 
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grange ■ -7.0 n. mi. 


( 6 . 22 ) 


The amount of range increase over a pure Breguet cruise as a result of 
optimal transition from cruise to a lower altitude is 


^Ib 


(6.23) 


A range » +7. 6 n. mi. 


(6.24) 


The sum of the increments is +0.6 n. mi. as has already been shown. 


The transitions will of course require less fuel than the 

(C^;.)^„„ - (C^)^. transitions. The amount of this fuel saving is calcu- 
r* mdx r min 

lated from c^ and c^.. In ascent 


AW » Ac.W.c « 466 lb. 
and in descent 

* y 1 

AW * ACfW.c e . 54 lb. 


(6.25) 


(6.26) 


for a combined weight saving of 520 lbs. The minimum value of Cp was 
taken as 0.011 instead of 0.012 so that the entire flight could be made 
with the afterburner on, that is, with a uniform engine description through- 
out the flight. 


Figure 6.8 shows the additional range realized by using u^p^ instead 

of (Cr) in an initial transition to cruise. The increment in range is 
r max 

plotted as a function of the initial value of with the initial value of 

A 

AHq taken as -0. 30. Incremental savings in mass, from Fig. 6. 1, are 
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converted to range savings by the modified Breguet range equation. The 
additional range in these transitions is comparable in size to the increase 
claimed (3. 8 n. mi. ) for complete trajectories with climb and descent 
transitions. in [ 10 ] . 

The demonstrated fuel savings are. of course, a small part of the 
total weight of the aircraft. Indeed, the use of singular perturbation 
methods implies that weight saving relative to the total weight will be on 
the order of c in comparison to one. So will the resulting increase in 
range when compared to the total range. But as a percent of payload the 
saving is not insignificant, since the percentage of payload to gross weight 
for an SST may be only on the order of 5% [ 25 ]. Furthermore, flight 
experience with the first operational SST, the Concorde, has shown it to 
have fuel reserves only on the order of 24, 000 lbs after a flight of 3400 
n. mi. (equivalent to a Paris to Washington, D. C. flight) carrying a 
payload that also happened to be 24, 000 lbs [ 26 ] . 

6. 4 McDonnell Douglas F- 4 

We next consider an early version of the McDonnell F-4 (Appendix 
A). This aircraft is capable of supersonic cruise, but we shall consider 
it only on transonic cruise. The principal reference for this aircraft (1) 
assumes that it has constnat SFC. The authors recognize a weak- 
ness in their assumption but justify it on the fact that better data were not 
available to them. We shall use this aircraft to observe the effect on cost 
of various values of the parameter y which is proportional to the slope of 
the curve of SFC vs Cp in the vicinity of the cruise value of SFC. We 
assume that the cruise value of SFC does not change as y changes. There* 
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fore, since 


SFC Cq + CjCp 


= c. 


1 + y 


D 


o J 


and to zero order in cruise 


we have 


SFC * Cq[ 1 + y(l + x^)l 


(6.27) 

(6.28) 

(6.29) 

(6.30) 


The value of x is determined solely from y. Since 
remain constant, the value of c^ must change with 


affect the value of e since 

c„ 


SFC at cruise is to 
y. Changing Cq will 


(6.31) 


pv 


Finally, from the equation for SFC, we see that the maximum value of 
SFC will increase as y increases. 

Figures 6. 9 through 6. 13 show a comparison of the cost between a 

full throttle climb to transonic cruise and the optimal cost for five values 

of y. First notice that if y is zero the cost of both trajectories is less 

than it would be for any other value of y. Then as y increases, the cost 

of both trajectories increases but the difference between them becomes 

greater. The largest value chosen for y is slightly larger than the value 

based on the Boeing SST data. The middle value of y corresponds to the 

cruise value of the Boeing 707- 320B which uses the PW JT3D turbofan 
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Figure 6. 9. Fuel Cost for u^^^ and Ini 

Various Initial Conditions at y » I 
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Figure 6.10. Fuel Cost for u^^p^ and Initial Transitions for 

Various Initial Conditions at y - 0. 01 










engine. The effect of Increasing y Is to Increase the value of SFC at 
maximum thrust. As y Increases, the cost of operating above cruise 
‘.hrust is increased. The optimal trajectories therefore tend to become 
much more heavily damped as y Increases. This relationship of y to 
damping was mentioned in Chapter V. Figure 6. 14 is a comparison of 
the transient responses in ^\s a function of Rj^. It shows clearly the 
relationship between y and damping. 

Figure 6. 15 shows the fuel saving in pounds associated with t!ie 
various values of y. It is a restatement of the data of Figs. 6. 9 through 
6. 13 for an initial weight of 30,432 lbs. It is seen that if y is zero the 
weight saving is only on the order of 10 to 20 lbs. However, for larger 
(but not unrealistic) values of y substantial savings in fuel can be realized. 


In Fig. 6. 16 the effect of the parameter y on the range improvement 
for the u^pj initial transition over that for is shown. The fuel 

weight savings from Fig. 6. 15 are converted to increments of range by 
the Breguet relation 


grange 4»»(1 - A M)(-;^)(— ) n. ml. 

^ Cq B07b 

^Ib ” 

The first order cruise flight path angle has been shown to be 


^b ’ 



(1 x^) 

(3x^ - 1) 


(6.32) 


(6.33) 


The parameter c^ changes with y so that SFC at cruise is constant. Its 
value is also plotted in Fig. 6. 16. Fuel savings are based on the data of 
Fig. 6. 15 at an initial flight path angle of zero. It is seen from Fig. 6. 16 
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that for values of y that are almost all less than that of the SST (i. e. , 
0.357), the F>4 achieves range increases that are generally much better 
than that of the SST. 

The result of this parametric study is to enable us to describe the 
conditions under which maximum range transitions to cruise can produce a 
worth-while saving in fuel. First it has been seen that the thrust re- 
quired in cruise must not be too close to the marimum thrust capability 
of the aircraft. If required cruise thrust approaches maximum thrust, 
fuel saving will approach zero. Second, the specific fuel consumption 
associated with maximum thrust must be greater than that required for 
cruise. If it is not, then the potential fuel saving will be negligibly 
small even though the maximum thrust may be much greater than the 
required cruise thrust. In summary there must be sufficient thrust 
capability over and above that required for cruise and there must be a 
cost associated with using it. 

A large number of aircraft do meet these conditions but it is also 
important to note that a large number fail to meet these conditions. Most 
of the commercial aircraft currently in service with the airlines fly at 
nearly their maximum thrust and in a fairly flat part of the curve of SFC 
vs Cp. This holds for the PW JT3D which powers the Boeing 707 and 
the McDonnell Douglas DC- 8, and also for the PW JT8D which powers 
the Boeing 727 and 737 and the McDonnell Douglas DC- 9 [ 29 ] . Both of 
these engines have a value of y on the order of 0. 04 and cruise at approx- 
imately 80% of maximum thrust. For the 707 in an initial transition to 
cruise, there appears to be no first order difference in cost between the 

Uopt policy and any of the other policies. For example, from an initial 
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State of (-.30, +. 10) space, the fuel cost associated with optimal 

transition to Breguet cruise is 0. 360cW.. The costs of the other three 

policies are 0.362cW., 0.372CW.. and 0.374cW. for (C .>) . , (C..) 

^ i* i* i r cruise' r max* 

^const* respectively. 

6. 5 Implementation of the Optimal Control Policy 

This thesis has developed the optimal control policy in terms of 
incremental changes in lift away from its cruise value. The control thus 
developed is a linear combination of the elements of the state vector and 
hence is suitable for implementation as a feedback control. It is also 
possible to implement an open- loop or programmed control system. Since 
analytic solutions are obtainable it is necessary only to program one of 
the zero- order state varibales, or AHq, or the control variable, u^, 
as a function of range. It would even be convenient to program which 

is as a function of range. 

It is probable, however, that to follow closely an optimal transition 

A 

trajectory would require a degree of accuracy in the measurement of AKq, 

AHq or Yq that could only be achieved by an inertial unit. Certainly the 
transition described for the SST in Fig. 6. 5, requiring an overshoot in 
altitude of 624 ft in a transition requiring 60 n. mi. and two minutes to com- 
plete, could probably not be duplicated by the pilot using a clock and a rate 
of climb meter for his cues. Exact duplication, however, may not be required. 

Simulation of transitions with a pilot in the loop would be necessary 
to determine how well a pilot could follow an optimal trajectory, what cues 
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he would require, and how sensitive the cost Is to deviations from the 
optimal. In the case of the SST, pilots have encountered difficulties in 
trying simultaneously to maintain constant Mach number and to level off 
to cruise altitude [24 ]. These difficulties were observed both in ground- 
based simulations and in flight simulations of the SST climb profile. Pilots 
had difficulty in avoiding overshoot in altitude and in Mach number while 
monitoring Mach error and pitch attitude. No data were available on fuel 
cost associated with the overshoots. 

Cost savings achievable thro'^gh optimal transitions would have to 
be weighed against the dollar and weight cost associated with the cues 
needed to implement the optimal transitions. For an aircraft that already 
has an inertial navigation system on board, it would be a simple matter to 
implement the optimal cruise transition policies of this thesis. Individual 
cost determinations would have to be made for other aircraft. 



CHAPTER VII 


CONTRIBUTIONS, CONCLUSIONS. RECOMMENDATIONS 
FOR FURTHER WORK 


7 . 1 Contributions 

This thesis has contributed an approximate analytic solution to one 
of a class of airplane performance optimization problems for which even 
computational solutions have been extremely difficult to obtain [8, il j . 
The analytic solution was obtained by the use of singular perturbation 
techniques in conjunction with the minimum principle of Pontryagin. 
Solutions were obtained in cruise and in transitions to and from cruise 
and then these three distinct segments were matched asymptotically. 

Inclusion of the normal acceleration equation made possible the 
inclusion of maneuvering lift effects on induced drag, an effect prcvimisly 
appearing only in computational solutions. Singular perturbation tech- 
niques allowed the drag force, which can not be linear iz 3d in any meaning- 
ful way, to be expressed as a quadratic function of state and control vector 
elements. The costate for mass, which is associated with these quadratic 
terms in the variational Hamiltonian, was shown to be a constant, to zero 
order in e, thus producing a solvable linear- quadratic optimal control 
problem. 

Optimal control laws have been developed for constant velocity 
transition to and from cruise in three dimensional state space (altitude, 
flight path angle, and mass), and expressions for the cost associated with 


them have also be^n developed. These cost expressions can serve as a 
lower bound for purposes of evaluating other techniques of performing 
these transitions. 

7 . 2 Conclusions 

Two basic conclusions can be drawn from this thesis, one from the 
analytic point of view and the other from the practical point of view. 

The first conclusion is that singular perturbation methods offer the 
possibility of approximate analytic solutions to certain optimal control 
problems that would otherwise have to be solved by computational methods. 
An important class of such problems is airplane performance problems in 
which equations involving aerodynamic drag (mass and velocity equations) 
can be neglected as a zero order approximation to the solution in the 
neighborhood of the singularities. The analytic results should be useful in 
themselves but should also provide clues toward finding methods of easing 
the computational difficulties associated with higher order versions of 
these problems. 

Second, for some aircraft the optimal control policy for transition 
developed in this thosis offers the possibility of significant fuel savings. 
These are aJrcraft that do not cruise at or near their maximum power 
setting (including transonic cruise for aircraft that are capable of super- 
sonic cruise) and for which the specific fuel consumption at maximum 
p ■ .er setting is somewhat greater than that for cruise. 



7. 3 Recommendations for Further Worl 


The first recommendation would be to apply the techniques demon- 
strated in this thesis to the solution of other similarly structured optimal 
control problems, for example, problems in which velocity varies slowly 
in the boundary layer and is constant to zero order in c. 

The solution obtained in this thesis should be stodled in conjunction 
with higher order computational solutions to try to gain insights into the 
nature of the optimal solutions and thereby to determine how best to 
approach computational solutions to this and similar problems with higher 
order state vectors. 

One could also expand the present study to consider the nature of the 
control if the slope of SFC as a function of Cj, were not merely a straight 
line in the vicinity of cruise but a series of connected straight line seg- 
ments. Such a representation would be mors accurate in the case of most 
power plants. The value assigned to the slope in a linear approximation 
directly affects the damping Inherent in the optimal control through the 
parameter y (Fig. 5. 2). and y determines x. the altitude parameter which 
represents the difference between optimal cruise altitude and the altitude 
for maximum lift- drag ratio. Since the altitude for maximum lift- drag 
ratio remains constant, a change in x represents a changes in the optimal 
cruise altitude. Thus as d(SFC)/d(Cp) changes discretely from one value 
to another, so also do two Important aspects of the optimal control prob- 
lem: damping in transition and optimal cruise altitude. There appears 
to be no point In using higher order functions of Cp to describe SFC. If 
a quadratic function were used, then the mass equation would be third 



order In thrust. The third order effects would be lost In forming the 
linear- quadratic optimal control problem. 

The descent from cruise could also be studied further. The poss- 
ibility of decelerating flight at minimum thr\:st or at zero thrust has not 
been considered here. Nevertheless the cruise- glide solution Is a very 
real possibility for maximum range flight. Optimal gliding night Is 
easily described by energy state methods. The transition from optimal 
cruise conditions when the engines are shut off to optimal gliding flight 
could possibly be set up as a boundary layer problem. Higher dimensional 
glides might be established as perturbations about the single variable 
optimal glide from energy-state methods. 

It Is also possible to study the maximum range problem from take- 
off through cruise and to landing by combining energy-state methods for 
acceleration and climb and for deceleration and descent with the optimal 
transitions and cruise developed In this thesis. This would require 
patching of solutions as opposed to matching. The energy climb (computa- 
tional solutions) would be followed until cruise velocity Is attained. This 
state would become the Initial condition for transitions to cruise. Similar- 
ly the state at the end of cruise when the engines are shut down becomes 
the Initial condition for the transition to optimal glide. 

The nature of the most general optimal transition from cruise to 
descent would In Itself be an Interesting study. Is there a throttling 
solution that Is superior to an Instantaneous zeroing of thrust? Should 
transition from cruise to the higher altitude for optimal glide be made 
while throttling or at constant thrust (Including zero)? 


It it hoped that this thesis will stimulate others to pursue these 
and other related topics of research In optimal airplane performance 
and the application of singular perturbation techniques thereto. 



APPENDIX A 


CHARACTERISTICS OF TYPICAL SUPERSONIC AND 
TRANSONIC CRUISING AIRPLANES 

This appendix presents those parameters necessary to describe 
the airframe and power-plant of three aircraft that cruise at constant 
Mach number in the stratosphere. 

A. 1 Boeing SST 

The first is the Boeing supersonic transport. Data is from [ 10 ]. 
For the airframe we have 


'^l 

640, 640 lbs 

initial cruise weight 

S 

7578 sq ft 

wing area 

V 

2479 ft/sec 

cruise speed 

M 

2.56 

cruise Mach number 


0.00878 

profile drag coefficient 

K 

0.5 

induced drag coefficient 


Table A. 1: Boeing SST Airframe Characteristics 

The power plant consists of four turbojets with afterburners, and in 
supersonic cruise the afterburners are on. Specific Fuel Consumption 
as a function of thrust coefficient is shown in Fig. A. 1. As a linear 
approximation to this function in the afterburning region we can write 

SFC • Cq + CjCp (A, 1) 
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SFC ■ 0.000265 + 0.010789 Cp 


(A. 2) 


with dimensions of Inverse seconds. 


Now certain parameters which are defined In Chapter III can be 
evaluated for the SST. If we write 


SFC » 1 + y 


"O ° 


then we have 


y - 0.357 


Since X Is determined solely by y we have 


x(y) ■ 0.762 


The first order value of the cruise flight path angle becomes 


Cjj . K) - 0.2151 


The cruise value of Cp becomes 


Cp(x, Cp ) - 0.01388 


which Is 60% of Us maximum value. 


153 



Finally, the parameter t becomes 


<(v, Cq. 3) - (1/450) (A. 9) 

In the stratosphere the scale height of the atmosphere is 

* 20,800 ft (A. 10) 

The range-optimal control for this airplane is a damped sinusoid 
in R. It has the following natural frequency and damping ratio, 
respectively 

= 0.33 (A. 11) 

Z = 0.8762 (A. 12) 

The parameter a, which appears in the y equation, is equal to the square 
of 

a = 0. 1089 (A. 13) 

A. 2 McDonnell Douglas F- 4 

The second airplane is an early version of the McDonnell Douglas 
F-4. It was used by Bryson, et al [ 1 ] and recurs frequently in later 
literature. This aircraft is capable of supersonic cruise but we consider 
it in transonic cruise. For the ainrame 
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^i 

30,452 lbs 

Initial cruise weight 

s 

530 sq ft 

wing area 

V 

872 ft/sec 

cruise velocity 

M 

0.9 

cruise Mach number 


0.014 

profile drag coefficient 

K 

0.2095 

induced drag coefficient 


Table A, 2; F-4 Airframe Characteristics 


The power plant consists of two GE J-79 turbojets with afterburners. 
Specific Fuel Consumption is taken nominally as 0.000625 sec”^. That is, 
if 


SFC 



(A. 14) 


then y is assumed to be zero. In the parametric study of Chapter VI^ y is 
allowed to assume various constant values while SFC at cruise remains 
constant. Figure A. 2 shows the extremes of this function. It is obvious 
that Cq, the Intercept on the SFC- axis, changes and hence the perturba- 
tion parameter, c, changes too. The cruise value of Cp is not allowed to 
vary. It holds constant at 0. 028, which is 32% of its maximum value. 


The parameters of Chapter III all depend on the value assigned to y. 
As a result they are presented in Table A. 3 for various values of y. The 
roots of the range- optimal control are shown in Fig. A. 3. 
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Figure A. 3. Foots of u . Equation for 


-• I Iipiw| 


y 

0 

0.01 

0.0425 

0.2 

0.4 

X 

1 

0. 9813 

0.9336 

0.8165 

0. 7517 

yib 

0. 1083 

0. 1104 

0.1172 

0. 1474 

0. 1833 

ISOOCq 

1 

0. 9807 

0.9263 

0. 7500 

0.6150 


67.0 

68.0 

72.0 

89.0 

109.0 


Table A. 3: F-4 Cruise Parameters for Values of y 
A. 3 Boeing 707-32QB 

The third airplane is the intercontinental version of the Boeing 707. 
The airframe parameters are from [27 ] , with drag coefficients extracted 
from Information on cruise thruse in [ 27 ] and [28 ]. We have 


W. 

i 

270,000 lbs 

Initial cruise weight 

S 

2892 sq ft 

wing area 

V 

775 ft/sec 

cruisa velocity 

M 

0.8 

cruise Mach number 


0.0114 

profile drag coefficient 

K 

0.062 

induced drag coefficient 


Table A. 4: 707-320B Airframe Characteristics 

The power plant consists of four PW JT3D turbofan engines without 

duct-burning capability. Specific Fuel Consiunption as a function of thrust 

coefficient is shown, for cruise in the stratosphere, in Fig. A. 4. The 

plot was developed from data in [ 20 ] . The aircraft cruises in the 
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SFC 



0.15 


aio 


h = 35,000ft 
and h > 40,000ft 


0.05 


0 


I L 

0 0.01 


.J-. 

0.02 


J I ( 

C.03 0.04 


Figure A. 4. SFC vs Cp for JTS-D 
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positive slope region close to maximum thrust coefficient. Note that a 
range- optimal solution will not use the negative slope region which would 
Increase SFC as Cp is decreased. Instead the range optimal solution 
would resort to chattering: minimum SFC would be maintained as thrust 
was reduced below the value for SFC^^^^ by alternately using zero thrust 
and thrust for SFC^.^. The duty cycle would be determined by the amount 
of thrust required. 

As a linear approximation to the function of Cp In the vicinity of Its 
cruise value we can write for SFC, from (A. 1) 


SFC * 0.0002014 + 0.0007508 Cp (A. 15) 

with dimensions of inverse seconds. The parameter y is 

ylCo* Cj, Cp ) » 0.0425 (A. 16) 

o 

The other parameters of Chapter III are 


x(y) = 0. 934 


(A. 17) 


y jj(x. Cp , K) = 0.0575 
o 

and 


(A. 18) 


«=F>cruis. " 

which is 80% of Its maximum value. Finally, the perturbation parameter 
Is 


€ » (1/185) 


(A. 20) 


60 


The range- optimal control for this airplane has two real roots 


f 


The values are 


Pj » -7.493 (A. 21) 

Pg = -0.149 (A. 22) 


The parameter a has the value 

a ■ 1.116 (A. 23) 


i 


? 


I 

I 

i 

i 
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